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ABSTRACT 


The  enclosed  development  study  was  perfocrr.ed  on  the  Bureau  of  Naval 
Weapons  Contract  NOw  61-0355-d. 

Phase  l  of  this  study  was  initiated  to  develop  a  "Generalized  Honeycomb 
Thermal  Stress.  Analysis'  which  would. after  the  capability  of  predict 
thermal  limitrfor  any  aircraft  subjected  to  the  effects  of  a  nuclear  explosion. 
The  method  developed  for.  this  contract  is  based  upon  a  recently  developed 
"AfTowabfe  Stress  Method*  with  the  capability  of  predicting  the  allowable 
load  on  any  structural  cross-section  in  any  flight  and  thermo!  exposure 
environment.  The  "Allowable  Stress  Method"  is  considered  one  of  the  most 
important  developments  in  analysis  of  aircraft  structures  in  recent  years. 

Such  an  approach  on  this  contract  gives  the  user  not  only  a  nuclear  effects 
analysis  method  but  a  "Generalized  Allowable  Structural  Loads  Method". 

Phase  II  of  this  study  was  initiated  to  develop  design  procedures  for  a  wea¬ 
pons  system  subjected  to  the  effects  of  a  nuclear  explosion.  The  design 
procedures  developed  during  this  study  provide  a  method  of  deriving  the 
minimum  weight  configuration  to  provide  the  needed  delivery  or  escape  capa¬ 
bility.  The  optimization  of  composite  structures  with  a  thermal  gradient 
presen*  is  a  break-thrcugh  in  automated  design  methods.  The  procedure  with 
minor  modification  can  be  used  for  the  aerodynamic  heating  problem  as  well 
as  the  nuclear  effects  heating  problem. 
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Specific  heat  of  core,  CaS/fpm-^f 
Specific  heat  of  face,  Cal/gmr^f 
Young's  modulus  for  the  core,  fc/in? 
lower  cop  strain,  m/m» 

Modulus  of  elasticity  at  temperature  of  element 

V%  Ibv'in2 
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Reduced  modulus  =  2E  Ej/Et  +  E 

Secant  modulus  =  f/i  +  .428  {F^,FCy)*t_^ 

Tangent  modulus  =  i/1  +  .428  m  (Fc/Fcy),n"’^ 

Elastic  applied  load  strain,  ii/tn 

Critical  column  buckling  strain,  It/in; 
critical  buckling  strain  of  plate  element,  ir/in 

Lower  panel  strain,  ir/in 

Bending  strain  on  lower  panel,  i**/in 

Total  strain  in  element  "n",  ii/in 

Element  strain  denoting  new  origin  of  stress-strain 
curve,  ir/in 

Axial  strain  for  inelastic  effects,  ir/in 
Strain  in  element  "nH  nt  temperature  step,  ir/in 
Spar  web  strain,  ir/in 
Applied  elastic  stress,  llM/in2 
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Bending  stress, 

A 

Compression  stress,  Ibs/tn 

Allowable  crippling  stress, 

Crf*-off  stress  for.  crippling,  Iby^in^ 

Stress  in  element  "n",  Ibt/h?' 

Ultimate  shear  stress,  Ibs/in? 

Ultimate  allowable  tensile  stress,  lbs/in^ 

*  Tension  or  compression  yield  stress  of 
element  "n“,  Ibs/in^ 

Yield  stress,  Ux/in^ 

sy 

Shear  modulus  of  core  material;  lbs/inz 
heat  capacity  of  honeycomb  core,  Cal/c it?  °F 

Heat  capacity  of  facing,  Cal/cm^  °F 

Heat  transfer  coefficient,  Cal/cm^  sec  °F 

Upper  stringer  deptb,  In 

Moment  of  inertia,  in* 

Total  elastic  moment  strain,  in/in 

Elastic  applied  moment  strain  about  the  x-axis,  in/in 

Compression  buckling  coefficient;  thermal  conductivity, 
Cal/cm  sec  °F 

Compression  buckling  coefficient  for  element  "n" 
depending  on  edge  fixity 

Moment  strain  term  for  initial  deflection,  \n/tn 

Internal  moment  strain  about  x-axis,  i n/in 
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SYMBOLS  (Continued) 

CefTsfze,  fn 
Temperature,  °F 
Core  foil  thickness,  is.  . 

Facing  thickness,  In. 


Thickness  of  element-  *Vk,  frr. 


Lower  stringer  thickness,  in. 

Upper  skin  thickness,  fn. 

Web  thickness,  In. 

Lower  skin  thickness,  fn. 

Lower  stringer  thickness.  In. 

Base  temperature  /  °F 

Deflection  of  structure  from  its  original  unloaded  position, 
in..  Spar  spacing.  In. 

Maximum  deflection  at  center  of  plate  element,  in. 

Distance  from  vertical  elastic  neutral  axis  to 
centroid  of  element  V  t  'n 


Distance  from  vertical  reference  axis  to  centroid 
of  element  "n",  in. 

Distance  from  horizontal  elastic  neutral  axis  to 
centroid  of  element  "n",  in. 

Distance  from  horizontal  reference  axis  to  centroid 
of  element  "n",  in. 

Thermal  coefficient  cf  expansion  of  element  "n",  in./in./°F 
Radius  of  gyration,  in. 
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SYMBOLS  (Continual 


Poisson's  ratio. 


Incremental  change 

Honeycomb  thermal  -response  parameter  • 

Density  of  core  1x,/in.^ 

Density  of  facing-  Ifer./fn^ ' 

Thermal  gradient  on  honeycomb  panel  °F 
Temperature  rise  on  front  face  of  honeycomb  panel  °f 


Time  to  peak  irrodl  once,  sec. 

-  K+fe  *  Vy)F 

.  1  +  L+VxVy  F 

TT- 


=  C1  +  2C2  +  C3 

♦  (q+JjJL-C,)  ♦  (c3+I^Lcjvy 

-  c,c3-cf  *  l^L 

-  a  tpj  tF2  E1E2 

(tpj  El  +  *F2E2)  Gcx  0“  ^ 

“  2  C  >Fi  tF2  eI£2 _ 

Opi  e1  +  fF2E2)  °2  Gcy  0  “  h 

=  effective  shear  modulus  of  core  in  x  direction#  Ih^oJ2 
=  effective  shear  modulus  of  core  in  y  direction#  fb^in? 
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SYMBOLS'  (Continued!} 

TKe  constants  Cj,  C^,  C%,  Cj  ore  dependent  upon  the  edge  fixity  < 
dl t ions  of  the  panel.  A  fisting  of  these  values  follows: 

Edges  simply  supported. 

Cj  *  CA  m  ^  —  1  C  j  *  *?(? 


n?o^ 


h  —  panel  length#  *» 
n  =  number  of  haff  waves 


*= panel  width,  hr 


Loaded  edge  simply  supported 
Unloaded  edge  clamped 

Ci  =  5.33  b2  C2  =  1.33  C, 

1  n2aZ"  1  3  b2 


C4  =  1.33 


7?" 


Loaded  edge  clamped 
Unloaded  edge  simply  supported 

c2=  i  c3=r'W.+  i).  i* 

^?  +  1  b2 

For  n  =  1  C]  =  .75  t2  C4  =  3  b2 


For  n  >  1 


Cl  =  C4  “  C-y-rr)  £ 


All  edges  clamped 

C2  =  1 .33  C3  =  — 

*  JTTl  b2 

For  n  =  1  Ci  =  4C  #  =  4 


Forn  >1  C.  =  4CA  =  5.33  b2 

(n2  +  1)  a2 
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SUBSCRIPTS 

cr 

AppfTerf  JofaT 

°p 

Applied 

cc 

Crippling 

c.o. 

Col-off 

cr 

Critical 

CO 

Compression  ultimate 

«y 

Compression  yield 

« 

• 

Effective;  equivalent 

i 

Initial;  internal 

• 

i 

Step  in  temperature-load  cycling 

ID 

Middle;  extreme  fiber 

n 

Element  number 

P 

Inelastic 

P* 

Permanent  set 

r 

Iteration  step 

T 

Temperature 

X 

About  x-axis 

y 

About  y-axis;  yield 

M 

Column  element  number 

NORTH  AMERICAN  AVIATION  ,  INC. 

COiuMtut  imMtt 
uuwwi  w— m 


REFERENCES 

Thermal  Stresses,  by  8.E.  Gatewood,  McGrnwHttt,  1952- 

Aircraft  Structures,  by  D.J.  Peery,  Me  Graw-Hill,  1950. 

NA59-1957,  "Structural-  Development  of  Brazed  Honeycomb 
Sandwich  Construction  for  the  Model  F1Q8  Air  Vehicle". 

NA58H-446,  “Criteria  and  Methods  for  the  Analysis  of  the 
Nuclear  Weapons  Delivery  Capability  of  the  A3J-1  Weapons 
System"! 

MIL-HDBK-5,  "Strength  of  Metal  Aircraft  Elements". 

NA56H-525,  "Structural  Tests  -  Complete  Airplane  for 
the  Model  A3J-1  Airplane". 

Forest  Products  Report  Number  1583— B,  "Effects  of  Shear 
Deformation  in  the  Core  of  a  Flat  Rectangular  Sandwich  Panel". 

NAoOH-288,  "Proposal  for  a  Generalized  Honeycomb  Thermal 
Stress  Analysis". 

NA60H-633,  "Correlation  of  Calculated  and  Measured  Thermal 
Stresses  in  Representative  FJ-4B  Aircraft  Aluminum  Structure 
Exposed  to  Transient  Heating". 


NORTH  AMERICAN  AVIATION  ,  INC 

cotwatut  >i»im» 

COtUMtHt*  s*„  «M» 


1.0  INTRODUCTION 

The  "Generalized  Honeycomb  Thermal  Stress  Study",  was  Initiated  in  two  phases  4s 
establish  the  capability  for  analysis  and  design  of  aircraft  subjected  to  nuclear 
explosions.  Phase  I  of  this  study,  "Analysis  and  Tests",  is  discussed  in  Sections  2.0^ 
3.0,  4.0  and  6.0.  Phase  II,  "Optimum  Design  Procedures",  is  discussed  in  Section 
5*<L 

The  analysis  procedures  developed  during  the  Phase  I  studies  are  presented  in  two 
sections:  a  detailed  exact  method  and  an  approximate  hand  calculation  method. 

The  exact  method  Is  based  upon  a  newly  developed  "Allowable  Stress"  method. 

This  method  provides  a  tool  by  which  the  allowable  lead  for  o  given  structure  can 
be  computed  for  any  of  three  failure  criteria;  initial  instability,  yield  and  ultimata. 
The  initial  stability  failure  is  based  upon  any  critically  instable  element  in  tho 
structural  cross-section.  The  yield  and  ultimate  failure  modes  ere  conversely 
based  upon  the  yield  or  ultimate  of  the  total  section  and  not  of  any  one  given 
element  as  has  been  the  case  In  the  past.  The  load  carrying  capability  of  each 
element  during  post-failure  as  well  as  pre-failure  loading  is  accounted  far  to 
best  predict  the  yield  and  ultimate  of  the  total  structural  section.  This  exact 
and  complex  method  has  been  programmed  for  the  IBM  709  Digital  Computer  and 
demonstrated  on  a  typical  horizontal  stabilizer  for  correlation  of  the  theoretical 
method  against  test  results. 

The  hand  calculation  method  was  developed  from  the  normal  thermal  stress 
analysis  approach  with  the  elimination  of  as  much  of  the  detail  as  possible  without  * 
completely  destroying  the  accuracy  of  the  method.  This  method  in  more  detail 
was  used  in  the  past  to  predict  the  thermal  limits  for  the  FJ4  and  A3J  airplanes. 
Limits  computed  from  this  abbreviated  method  can  be  expected  to  be  as  much  as 
+  15%  in  error  depending  upon  the  core  token  while  making  the  analysis.  Froni- 
the  stress  distribution  computed  with  the  thermal  stress  equations,  a  comparison 
with  element  allowables  must  be  made  to  determine  the  limits  os  demonstrated 
in  the  example  problem  presented  in  Section  3.2.  Although  this  procedure  can 
be  accomplished  by  hand,  it  is  particularly  adaptable  to  desk  calculators  and 
small  digital  computers. 

The  Phase  I  study  also  Included  high  temperature  structural  tests  of  a  typical 
horizontal  stabilizer  and  several  honeycomb  panel  configurations.  The  typical 
horizontal  stabilizer  test  was  made  to  determine  an  ultimate  failure  load  under 
a  transient  temperature  condition.  This  failure  load  was  then  compared  with  a 
computed  load  showing  excellent  agreement. 

The  honeycomb  panel  tests  performed  during  this  study  established  test  data 
for  correlation  with  a  general  honeycomb  panel  buckling  equation.  Previous 
to  this  study,  honeycomb  panel  buckling  allowables  were  determined  by 
testing  only.  With  the  correlation  of  the  generalized  buckling  equation. 
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allowables  for  all  the  structural  elements  can  be  computed  theoretically, 
thus  prov uling  the  capability  of  analysis  and  design  far  thermal  effects ae 
structures  without  excessive  testing. 

The  Phase  II  study  has  approached  the  objective  of  aircraft  design  for 
induced  thermal  gradients  os  a  design  optimization  problem.  In  general, 
the  nuclear  effects  problem  has  been  handled  by  determining  the  limit  which 
the  aircraft  can  withstand  after  the  design  has  been  completed.  For  the  oqm 
of  aircraft  vulnerability  to  defensive  weapons,  a  capability  must  be  designed 
into  the  -ilrcraft.  In  many  cases,  it  would  also  be  a  great  advantage  to  pre¬ 
determine  the  delivery  capabilities  os  the  aircraft  is  designed,  fnitia re¬ 
designing  for  the  required  capabilities  will  insure  the  best  design  for  the 
weapon  system. 

The  design  optimization  methods  developed  during  this  study  include;  1. 
honeycomb  panels  with  a  fixed  load,  2.  honeycomb  panels  with  a  fixed 
strain,  3.  honeycomb  box  beam  with  straight  or  corrugated  spar  webs,  4. 
skin-stringer  box  beam  with  straight  or  corrugated  webs.  In  each  of  these 
cases,  a  minimum  weight  design  is  computed  from  a  given  load  and  thermal 
input.  The  temperature  of  the  hot  facing  is  a  variable  dependent  upon  the 
materia!  thicknesses.  The  following  dependent  parameters  Ore  also  determined 
in  these  optimization  procedures;  honeycomb  facing  thicknesses,  core  depth, 
core  cell  size  and  foil  thickness,  spar  spacing,  spar  web  thickness,  stringer 
spacing.,  stringer  thickness  and  depth. 

The  methods  developed  in  this  study  provides  a  basic  tool  for  nuclear  effect 
capability  design.  They  have  also  established  a  basis  from  which  further 
design  methods  can  be  developed. 
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2.0  GENERAL  DISCUSSION  OF  THE  ALLOWABLE  STRESS  METHOD 

When  a  structure  fs  subjected  to  a  non-uniform  temperature  distribution, 
the  elements  of  die  cross-section  of  the  structure  hove  varying  material 
properties  depending  upon  the  temperature  of  the  element-.  Thermal  stress** 
may  be  present  -due  to  a  fton-uniform  temperature  distribution  and/or  differ¬ 
ent  materials  in  the  cross-section.  The  inelastic  effects  of  buckling  and 
stresses  beyond  the  proportional  limit  produce,  stress  variations  and  dishibu— . 
tions  within  the  cross-section.  The  strains  associated  with  the  stresses  vary 
due  to  thermal  expansion  applied  axial  loadand/or  bending  moment  strains 
and  stress  variations  produced  by  the  temperature  distribution. 


With  the  material  property  stress-strain  curve  of  each  element  in  the 
cross-section  represented  analytically  by  the  Rambsrg-Osgood  equation 
(Ref.  1)  then. 


where,  for  element  n,  is  the  yield  stress  which  depends  upon 


temperature,  previous  temperature  and  load  history  of  the  element;  Fn 
is  the  stress,  En  is  the  elastic  modulus,  and  en  is  the  element  strain  associated 
with  the  stress. 


The  following  sections  show  the  development  of  the  strain,  stress,  and 
equilibrium  equations  necessary  to  define  the  load-deformation  curve  for  any 
cross-section  in  any  temperature-load  environment. 


2 . 


2.1.1  BASIC  STRAIN  EQUATIONS 

Assuming  a  plane  cross-section  to  remain  plane  at  all  time*/  then. the  shape 
of  the  strain  distribution  curve  on  the  cross-section  will  olwa/s  be  the  same 
shape  as  the  thermal  strain  distribution  curve  except  for  possible  iototTow 
through  a  constant  angle.  This  relationship  holds  true  wftetherthe  material 
is  elastic  or  inelastic  and  regardless  of  the  presence  of  buckling.  An  applied 
c.xial  lead  produces  a  constant  strain  and  an  applied  Moment  produces  <r 
ro.'ation  through  a  constant  angle.  A  general  strain  equation  may  be  written 
for  any  element  n  in  the  cross-section  based  on  the  above  assumption.  The 
element  strain,designated  as  en/  consists  of  the  sum  of  the  elastic  thermal 
strain,  the  elastic  applied  axial  load  and  bending  moment  strains,  and  a 
correcting  linear  term  to  accouni  for  all  inelastic  effects'.  This  strain  may 
be  written  as  * 


;n  =  eTn  +  ?apn  +ef 


(2.1.1) 


where  eT  is  the  elastic  thermo!  strain  and  is  defined  as  follows: 

•n 

eTn  =  -  («T)n  +  eT  +  Kjx  jy„j  +  Kj^  (2.1.2) 


where 


=  1  <°U„EnAn 

**n 
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Hie  tewpewiaw  moment  terra  in  Eq;  (2. 7.3)  ore  defiriedby 


Mj  *  T  (  a  t)n  En  An  yn 

x  n 

•  •  * 

%,  *  5  (  ^ir "  ^Vr  “rr 


J*  (2.K4) 


and  the  elastic  bending  stiffness  parameters  about -orthogonal  axes  x  and  y 

are  defined  by 


(0«>  -  £e„  a„  y„ 

n 

^  At  Xn 

<£fxy>  *  5  ^  A*  xn  ^ 


(2.1.5) 


if  the  bending  stiffness  of  each  element  about  Its  own  elastic  neutral  axis  Is 
neglected. 

from  the  arbitrary  reference  axes  xg  and  yg  the  values  xR  and  yn  are 
defined  by 


xn»  "  n  ^n  At  *nR 


X  E  A 

n  n  n 

yn  =  yn„.  n  En  An  Yr-R 

X  E  A 
n  n  n 


(2.1.6) 


The  elostic  replied  axial  load  and  bending  moment  strain  term  in  Eq. 

(2.1.1)  i* 

£?»  -  %  +  K°Px  (-£)  +  S  (-r-) 


(2. 
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where  the  efostfc  appfTecf  axiaf  footf strain  Tt 


°P 


%  E„  A 


12. 


and  the  efostfc  applied  bending  moment  strains  on  the  extreme  fiber  elements 
are 


K^x  “  c 

Mx  (Ely)  -  My  <0W)  1 

(El,)  (Ely,  -  (EIXy)2 

\  <Elx  -  (EIXy, 

N  D 

opy 

_(EIX,  (Ely,  -  (EIXy)2 

In  Eq's.  (2. 1.8)  and  (2.1.9)  P,  M*,  and  M y  are  the  applied  axial  load, 
bending  moment  about  the  horizontal  elastic  neutral  axis,  and  bending  moment 
about  the  vertical  elastic  neutral  axis,  respectively,  for  any  given  cross- 
section. 

The  inelastic  correcting  strain  is  determined  by 

(?) 

where  the  values  of  the  inelastic  correcting  strains  ep,  and  Kp  are 
defined  by  the  iteration  procedures  of  Section  2.1. 7such  that  the  stresses 
Fn  corresponding  to  the  strains  en  satisfy  the  equilibrium  equations 

P  =  ^  FnAnCn 
*-?  FnAn^nVn 

My  =2#  FnAnCnXn  (2.1 


where  the  effective  area  factor  for  local  buckling  of  element  n  as  defined 
in  Section  2.1.6. 
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2*1.2  TEMPERATURE  AND  LOAD  STRAIN  EQUATIONS 

Since  the  applied  load  may  already  be  present  when  the  temperature 
Is  applied  to  the  structure,  a  non-uniform  temperature  distribution  may 
be  present  when  the  load  Is  applied,  the  temperature  may  cycle  with  a 
steady  load  acting,  or  Loth  load  and  temperature  may  wary  with  time*  It  h 
necessary  to  write  the  strain  equations  fa  account  for  sequence  application 
and  removal  of  food  and  temperature.  These  equations  must  then  allow  for 
the  strain,  stress,  and  temperature  histories  af  the  structural  elements* 

Various  factors  involved  are  the  true  permanent  set  strains  at  the  end  of  any 
given  sequence  of  temperature  and  load  cycling,  strain  accumulation  pro¬ 
duced  by  temperature  cycling,  variation  of  material  properties  as  the  temp¬ 
erature  varies  with  time,  and  the  recovery  of  material  properties  after  a 
given  exposure  time  at  temperature.  Under  sequence  application  and/or 
removal  of  load  and  tenperature^unloading  of  some  elements  from  an  inelas¬ 
tic  position  or»  the  stress-strain  curve  may  occur.  For  unloading,  the 
element  strain  is  assumed  to  follow  a  straight  line  of  slope  En  (elastic)  to 
zero  stress  and  then  follow  a  stress-strain  curve  with  a  new  origin  when 
loading  in  the  opposite  direction.  The  shift  in  origin  for  any  element  is 
defined  as: 

(2.1.12) 


(See  Figure  2.1.) 

If  reloading  takes  place  from  ony  point  on  the  straight  line  or  from  an 
elastic  position  on  the  new  stress-strain  curve,  then  the  straight  line  is 
followed  up  to  the  original  stress-strain  curve  (for  strain  hardening  materials) 
and  the  original  stress-strain  curve  followed  for  larger  strains.  If  reloading 
occurs  from  an  inelastic  point  on  the  new  stress-strain  curve,  then  the  element 
may  follow  a  new  stress-strain  curve  from  a  new  origin  as  shown  by  (6)  in 
Figure  2.1.  The  effect  that  actually  occurs  depends  upon  the  amount  of 
yielding  that  has  occurred  in  the  reverse  direction.  Until  more  test  data  is 
available  to  establish  definitive  criteria  for  reloading,  it  may  be  reasonably 
assumed  that  a  new  stress-strain  curve  con  arise  whenever  unloading  occurs 
from  any  point  on  the  stress-strain  curve  where: 

Fn  =  0.80 

1Vn' 


(2.1.13) 
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TKs  corresponds  Ta  a  shift  in  origin  of  the  sfress-sfrotrr  curve  defTnecf  by 


if  the  stress-strain  curye  is  represented  by  the  Romberg -Osgood  equation  with 
m  =  10.  (See  Figure  2.0  where  reloading  to  the  old  stress-strain  curve  at  the 
maximum  permissible  origin  shift  defined  by  Eq.  (2.1. 14)  is  illustrated  by  (JJ.  -  - 
ft  is  apparent,  then,  that  with  various  elements  in  the  cross-section  loading, 
unloading,  or  reloading  at  any  step  in  the  temperature-load  sequence,  that  the 
strains  and  stresses  must  depend  upon  the  previous  steps  In  the  sequence  to  define 
the  proper  stress-strain  curve  and  origin.  If  at  each  step  in  the  temperature- 
food  sequence  the  strain  en  is  assumed  to  start  from  a  new  origin  as  defined  by 
Eq.  (2.1. 12),  then  the  element  strain  Eq.  (2.1.1)  ot  any  step  {  can  be  written 
in  terms  of  the  previous  step  ]-1  and  the  load  or  temperature  application  or 
removal  at  step  }  os: 


enj  =  Fnj-1  -(«T)nj+  (®T  +  eap  +  V| 

£n,H 

+  (KTX+Kapx+Kp^  J^+(KTy+KCpy+Kpy)  . 

I  c  I 

*  Fn  .  , 

.  Aej.  (2.1.15) 

S-l  1 

where  bn,]-]  *s  elastic  modulus  at-  step  j-1 .  Eq.  (2. 1.15)  is  general 
since  it  includes  both  the  thermal  elastic  strain  of  Eq.  (2.1.2),  the  elastic 
applied  load  and/or  moment  strains  of  Eq.  (2.1 .7),  all  inelastic  effects  which 
are  represented  by  ep,  KpX,  and  Kpy,  and  the  stress  effects  of  the  previous  steps. 

Since  the  calculations  are  performed  for  temperature-load  sequences 
where  the  effects  of  temperature  and  load  are  not  applied  simultaneously, 
the  general  strain  equation  (2.1.15)  is  written  as  follows.  For  any  temper¬ 
ature  application  or  removal  step;  the  element  strain  is: 
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*1  *  *n,H  ♦  &«h, 

K.. , 

ri 


"‘H  -  («rj r  +  feT+  e  l  +  (Kj  +  Kp  }  + 

1— -  i  P|  *  p*i  c 


(Kt  +  Kp)  *n_ 

y  y  b 


(2.1.16) 


In  Eq.  (2.1.16)  the  temperature  Tn.  is  the  incremental  change  In  temp¬ 
erature  of  element  n  from  the  last  temperature  application  or  removal  step. 

The  temperature  terms  er,  and  Ky.  are  essentially  as  shown  in  Eq.  (2.1.3) 
and specifically  for  any  lemperatufe  application  or  removal  step  are  calculated 
as  follows: 


'Tj  ■  (*  (°Va). 


n  I 


KT  * 

*! 


-  ?(aVnVnl  <  *  EnAnW. 

n  I  "  J  n  |  n  I 

<fiA^  (SvW>,  -  (5E»A"v">i 


-  £<«lWWj<5  £/„*„/„) 


5i  - 

.  ‘S£A>&  (IV A  - 


(X  E  A  x  Y  ). 
n  "  nnn| 


(2.1.17) 


For  any  load  application  or  removal  step  |  the  element  strain  is 
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*n  ‘ 

S-i 

=  Fn;  i  +  (e  +  e  ).  +  (K  +  Kp  ).  + 

nH  op  P|  ap^  rx  1  “ 

^  %*■%>.  _% 


(2.1.18) 


wnere# 


w 

K  (S„EnV5j  -  “/j  R!»AnV",j 

ci  - 

.(S  EnVn)j  <KV5j  -  <KA„V»>? 

.  b  pVj  (PA^j  -  H=j  (IEAVn>j 

'  '5s4  lSEnVS>j  -  (fvvwj? 


(2.1.19) 


where  P-,  M^.,  and  My.  are  the  incremental  loads  and  bending  moments 
applied  or  removed  at  afiy  step  j. 

In  Equations  (2.1.15)  and  (2.1.18)  en.  is  the  strain  directly  associated 
with  the  stress  Fn.  either  as  an  elastic  value  or  as  a  volue  on  the  stress-strain 
curve  upon  reloading.  In  this  case  the  strain  associated  with  the  stress  is 
taken  as:  ,  , 


en:  +  V  WhCre  =  i  en0i 
•  »  •  i=k 


(2.1.20) 


where  en0.  is  determined  by  Equation  (2.1.12)  but  is  not  included  in  the  sum 
unless  it  also  satisfied  the  condition  of  Equation  (2.1.14).  k  is  the  step  at 
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which  the  last  sign  reversal  of  e^j.  occurred,  and  j-m  is  the  lost  step  of 
which  satisfied  Equation  (2.1*.  14). 


The  stress  equations  of  Section  2.1.5  and  the  logic  table  in  Table  2.1 
ore  required  for  the  sebtion  of  the  strains  and  stresses  of  any  step  f  m  the 
temperature-load  sequence. 
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2.1.3  SHIFT  IN  ORIGIN  OF  STRESS-STRAIN  CURVES 

Figure  2. 1  illustrates  the  various  stress-strain  relationships  which  can 
arise  during  the  seance  application  or  removal  of  temperature  and  load. 

Each  element  in  the  cross-section  is  considered  to  behove  according  to  one  • 
of  these  conditions  of  any*  step  irr  the  temperofare-lood  sequence.  Analytically 
these  conditions  ore  described  by  the  Ramberg-Osgood  Eq:  (2.1.22)  and  the 
logic  table  of  Table  2»1  which  allow  for  the  stress  and  strain  histories  of  the 
elements. 

Curve  (1)  represents  loading  on  the  original  stress-strain  curve  From 
the  0,0  origin.  Assuming  no  change  in  material  properties  due  to  tempera¬ 
ture  change,  (2)  represents  unloading  of  the  element  from  an  Inelastic  point 
on  (1)  on  an  elastic  line  which  has  the  some  slope  (E,,)  as  the  elastic  portion 
of  the  original  curve  (1).  Unloading  of  the  element  may  continue  until  the 
element  begins  to  load  In  the  opposite  direction.  In  this  case,  a  new  stress- 
strain  curve  (3)  is  followed  which  has  a  new  origin  defined  by  the  offset  strain 
eIlQ  where  en  is  found  by  Eq.  (2. 1 . 12).  The  elastic  unloading  lines  (4) 
show  unloading  from  Nro  different  points  on  the  stress-strain  curve(3).  One 
elastic  line  is  similar  to  elastic  unloading  line  (2)  and  takes  place  from  an 
inelastic  point  on  the  offset  stress-strain  curve  (3).  Continued  unloading  of 
the  element « lay  produce  reloading  of  the  element  in  the  same  direction  as 
the  original  loading.  This  reloading  is  assumed  to  take  place  along  a  new 
stress-strain  curve  (6)  which  has  a  different  origin  than  either  (1)  or  (3). 

The  other  unloading  line  (4)  begins  at  a  point  on  (3)  defined  by  Fn  =  0.80  *Yn 
which  is  the  assumed  Imitation  for  reloading  back  to  the  orib:r»al  stress-strain 
curve  (1).  This  defines  o  maximum  permissible  offset  shown  by  Eq.  (2. 1 . 14) 
for  the  Ramberg-Osgood  curve  with  m  =  10.  Continued  application  of  strain 
on  the  element  produces  reloading  on  an  extension  of  the  original  stress-strain 
curve  (1)  which  is  designated  by  (7). 
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T.T.4  STGIV CUNVENTTON 

The  use  of  the  equations  of  Sec.  2.1.1  end  2.1..2  requires  that 
definite  sign  conventions  be  established  for  compatibility  throughout  the 
temperature-load  sequencing.  The  following  table  shows  the  terms  requiring 
algebraic  sign  definition  end  the  proper  signs. 


EQUATION  TERM  ALGEBRAIC  SIGN 


T 

(+)  if  temperature  increases  from  datum 
(-)  if  temperature  decreases  from  datum 

xn,  Yn 

(+)  dimension  toward  the  extreme  fiber 

element  having  the  highest  positive 
value  of  (  aT).  In  the  case  of  symmetrical 
temperature  distribution  and  geometry 
any  convenient  sign  convention  may  be 
sefected  if  the  applied  moment  values 
have  consistent  signs. 

Vm 

(+)  tension 
(-)  compression 

c,  b 

(+) 

«ap 

(+)  tension 
(-)  compression 

(+)  if  the  applied  bending  moment  M^  puts 
compression  on  elements  having  -yn 
values. 

(-)  if  puts  compression  on  elements  having 
-tyn  values. 

{+)  if  the  applied  bending  moment  My  puts 

compression  on  elements  having  -xn  values. 
(-)  if  My  puts  compression  on  elements  having 
+xn  values. 


NORTH  AMERICAN  AVIATION  ,  INC. 

COlUMtuS  DIVISION 

COtcr«tt/$  r*.  gwy 


2. 1 .4  SIGN  CONVENTION  (continued) 


EQUATION  TEEM  .  ALGEBRAIC  SIGN  • 


Mk 

*v 


(+)  Tension. 

(-)  compression 

(+)  if  moment  produces  compression  on  elements 
having  -yn  values. 

(-)  if  moment  produces  compression  on  elements 
having  +yn  values. 

(+)  if  moment  produces  compression  on  elements 
having  -xn  values. 

(-)  if  moment  produces  compression  on  elements 
having  -b<n  values. 


en'em 


(v)  tension 
(-)  compression 
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2.1.5  STRESS  EQUATIONS 

Since  all  calculations  for  the  cross-section  are  performed  using  the 
strain  equation  (2.1.1)  it  is  necessary- to  obtain  the  stress  Fn  associated  «4tir 
the  strain  en  from  a  stress-strain  curve  Tor  each  element  corresponding  to  the 
temperature  of  the  element.  The  stress-strain  curve  for  each  element  in  the 
cross-section  is  represented  analytically  by  the  non-dimensional  Rambergr 
Osgood  equation  as  shown  in  Reference  (1). 


en  = 


+  ~  /Fn\ 

fe) 


(2.1.21) 


where  the  slope  of  the  elastic  portion  of  the  curve,  the  yield  stress,  and  a 
shape  factor  define  the  stress-strain  curve. 

For  the  stresses  at  any  step  {  in  the  temperature-load  sequence  corres¬ 
ponding  to  the  strains  of  Equation  (2.1.16)  or  (2.1.18)  the  Romberg- Osgood 
Equation  (2.1.21)  is  modified  to 


where 


(°"i +  ’-i 


|-m 

=  0  or  S  e_ 


_L  i  + 


W 


(2.1.22) 


rn.  =  0  or  1 


according  to  the  stress  and  strain  histories  of  the  elements  with  the  values  of 
qnj  and  rn.  defined  by  the  logic  table  of  Table  2.1.  In  the  modified  Ramberg- 
Osgood  Elation  (2.1 .22)  En-  is  the  elastic  modulus  of  element  n  at  step  j, 
Fyn.  is  the  yield  stress  of  element  n  at  step  j,  and  t„  is  the  shape  factor  which 
depends  primarily  upon  the  material. 
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2.1.6  BUCKLING  AND  CRIPPLING  LOADS 


For  stresses  below  the  fccaf  buckling  stresses  of  the  elements  of  the  struc¬ 
ture,  the  effective  area  coefficients  as  defined  by  Cn  in  the  iteration  equa¬ 
tions  of  Sec.  2.1.7,  are  token  as  unity.  However,  after  some-elements  buckle, 
the  stress  end  strain  distribution  on  the  elements'  change  due- to  the  change  Sr 
length  of  the  buckled  element  from  the  deflection.  Since  some  part  of  the 
buckled  element  (one  or  more  edges),  does,  not  deflept,  thks tnim  in  this,  under- 
fleeted  portion  of  the  element  may  be  token  os  the  reference  strain  for  the 
element  andfor  the  strain  to  be  associated  with  the  strains  of  the  other  elements. 


Now  the  buckling  strain  of  the  element 


n  can  be  approximated  by: 


(2.1.23) 


where  the  buckling  curve  is  taken  to  be  the  same  as  the  stress-strain  curve 
(see  Figure  6  -  5  in  reference  1).  The  temperature  of  the  element  is  assumed  to 
be  uniform.  In  most  cases,  this  will  be  approximately  true.  If  not,  an  average 
temperature  of  the  element  can  be  used.  Define  the  reference  strain  for  the 
element  after  buckling  by: 


(2.1.24) 


where  bne  is  an  effective  width,  Kn  depends  on  the  restraints  after  buckling  of 
the  element,  and  en  is  given  by  Eq.  2.1.24  with  local  bending  of  the  element 
neglected  but  with  ep  including  the  orea  change  effect. 


From  Eqs.  (2.1.23)  and  (2.1.24) 


(2.1.25) 


and 


X  F  A  C 

eop  -  n  n  n  n 


X  En  An 


Cn  *  h  en  -  ei 


cm 


(2.1.26) 
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The  maximum  value  of  e  Is; 

op 


* 

Fcc 

/E  *aP  J  =0.88 

/E  ecr 

/— 

at  _JF_  _ 

p 

Fey 

\  Fcy  /ax. 

Vfq,  \ 

1  Kcr 

rcy 

E  e  .  Eec,. 

=  0.98  >^-2-  . 

<y  qr 

F  ol 

E  e 

Eecr  > 

F 

^7 

i 

S' 

u. 

11 

0.87 


0.98 

(2.1.30) 


Equation  (2.1 .30)  indicates  that  if  the  buckling  strain  is  less  than  the  elastic 
yield  strain  FCy/E  then  the  maximum  crippling  stress  occurs  at  the  elastic  yield 
strain  while  if  the  buckling  strain  is  greater  than  the  elastic  yield  strain  the 
maximum  stress  occurs  ot  the  buckling  strain. 

i 
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2.1.7  ITERATION  EQUATIONS 

The  fTnaf  strains  and  stresses  erf  any  step  f  fn  the  temperature-food  sequence 
are  obtained  by  the  iteration  of  the  strain  equations  (2.1.16)  where  the 
inelastic  effects  represented  by  epj,  Kp^.  and  Kp  .  are  determined  by 
axial  load  and  bending  moment  equilibrtdm  on  the  cross  section.  At  any 
temperature  or  load  step  {  the  strains  are  calculated  by  Eq.  (2.1.16)  or 
Eq.  (2.1.18)  withep?,  Kp*./  aid  Kpjy  taken  as  zero  for  a  first  nppfmctnvrf-r<yr. 
These  strains  are  useo  to  obtain  the  element  stresses  from  Eq.  (2.T.22)  and  the 
logic  tcble  of  Table  2.1.  If  inelastic  effects  are  present,  these  stresses 
will  give  values  of  M*,  and  My  in  Eq.  (2.1*11)  different  from  the  valuer 
of  ?Pj,  fMxj,  and  used  in  Eq.  (2.1.34). 

The  inelastic  effects  for  successive  approximations  are  defined  by  the 
following  equilibrium  iteration  equations  at  any  step  j  in  the  temperature - 
:oad  sequence.  The  axial  load  and  bending  moment  equilibrium  iteration 
equations  are  shown  below  for  any  step  r  in  the  iteration. 

(S)  ' =  (epi)  '  (e°pi)  *  f  \ 

(V)  ,  '(Vj  r-l'W  r_,  +  fS 
(Kp/i)  r”  (K"i)  r-1  '  (Sj),.,  +  2  S, 


where  the  external  applied  loads  are  represented  by: 
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section  are  calculated  at  each  iteration  step  r  by  the  following  equations:  ■ 
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The  iteration  continues  until  the  successive  aoproximatioris-of  Eq.  (2.1.33) 
produce: 

ePj  =  (e PjJ  r  =(epi)  r  l-  Kpx\  =(KpXi)  ^  and 
Kp  _  =^Py*j  r  =  a  sPec'^et^  tolerance. 

A  practical  range  of  tolerances  has  been  established  ate 


^  -  ^Pjjr_ i=^0.00003  to +  0.00006 

N  r  '  (S)r-l  =i  °-°0003  ,0i°-‘ 

(K4r  ‘  (S)r-l  =1°-°0003  ,£>1°-° 


00006 


00006 
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2.1.8  COLUMN  PROCEDURES 

Foe  a  structural  cross-section  wfiTcfi  contains  column  elements,  it  is  fTrst 
necessary  to  define  the  load-deformation  curve  of  the  column  element  if  a 
satisfactory  representation  is  to  be  mode  as  described  in  Sec.  2.1.9.  The 
Ramberg-Osgood  stress-strain  curve  (2.1.21)  adequately  describes  many 
elements  in  the  cross-section  except  those  subject  to  either  column  or  local 
instability.  The  load-deformotion  chorocteristics  of  members  with  local 
instability  are  definedby  the  methods  of  Sec.  7.T.6  where  the  Romberg— 
Osgood  stress-strain  curve  and  an  effective  area  which  depends  upon  the 
critical  buckling  strain  of  the  element-  define  the  load-deformation  chanoc— 
teristics.  The  column,  however,  may  have  an  unsymmetrical  temperature 
distribution  through  the  cross-section  or  a  secondary  bending  moment  may  be 
acting  on  the  column  element  thus  producing  deflections  normal  to  the 
neutral  plane  of  the  column.  An  axial  load  acting  on  the  column  element 
wili  affect  these  deflections  by  producing  o  bending  moment  due  to  the 
deflections.  If  symmetrical  bending  only  is  considered,  tbe  net  bending 
moment  strain  consists  of  two  parts  ,  one  due  to  applied  secondary  bending 
moments  and  one  due  to  deflection. 

Kax  -  ^OpX  +  KWx  (2.1.36) 


The  deflection  moment  strain  may  be  described  in  terms  of  eap  and  the 
deflection  W  as  follows: 

In  many  problems  the  applied  secondary  bending  moment  may  be  related  to 
the  applied  axial  load  as 

MaPx  =  PP°P  (2.1.38) 

In  this  case  Eq.  (2.1.37)  may  also  be  written  , 


(2.1.39) 
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tf  z  is  .the  variable  along  the  length  of  the  structure,  then  W  depends  on 
z  and  the  end  conditions,  while  p  may  depend  on  z  If  the  applied  moment 
is  variable.  For  simple  beams.  Reference- (2)  ghreer 


d7W  =  _M_  =  JK_ 

dz2  El 


(2.1.40) 


The  net  rotational  strain  of  any  cross-section  through  the  column  may  be 
written 


Kx=Kix  +  K.»x  ♦  Kapi 


+  (Ni- 


,  +  KP 


+  *2  ^x. 


(2.1.41) 


for  the  temperature-load  sequence  +  T,  +  P  where  K-  represents  the  InltTaf 
deflection  or  eccentricity  of  the  load.  In  Eq.  (2.1.41)  Kpx  is  a  function  of 
z  depending  upon  the  amount  of  inelastic  action  at  the  various  cross-sections 
throughout  the  span  or  length  of  the  column.  Eq.  (2.1.40)  con  be  written 


(2.1.42) 


where,  with  K  as  a  table  of  values  for  selected  values  of  z,  then  Kwx  can 
be  obtained  from  Eq.  (2.1.42)  by  an  area-moment  numerical  integration  or 
by  a  double  summation. 


If  Kapx,  Kjx  and  K;x  are  constant  along  the  beam,  then  Kpx  will  be  con¬ 
stant  except  for  the  effect  of  W.  Assume  KPx  is  constant  along  the  span  so 
that  Eq.  (2.1.42)  can  be  integrated  to  give 


Kwx 

f  (*) 


cos  q  (1  -  )  -  cos  q 

- 1 -  (compression) 


cos  q 
cosh  q  (1 

cosh 


)  -  cosh  q 


q 


(tension) 


(2.1.43) 
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A* 

for  a  simply  supported  beaw  cotuow 
let  e_  P  *  , 


2 

q  = 


^  (f 

so  that  the  maximum  deflection  in  £q.(2. 1 .43)  at  z  =  ^  »* 

(KwJ  «/* i  +  *a  ♦  (KT).  •  ♦  ^  I 

x  max.  I  x  px  *x  {-] 

Eq.(2. 1.41)nKiy  be  written 
*  ^  — ! 


Px- 


Kv  =  1  +f  max. 


i-1 


(2.1.4^ 


(2.1.45) 


(2.1.46) 


^x^opx+(KT^_1+  S  Kpx+Kp)) 


(2.1.47) 


The  temperature-load  sequence  of  primary  importance  for  defining  the 
column  load-deformation  curve  is  apply  temperature  (+T),  then  apply  load 
(+f).  For  any  given  temperature  distribution  through  the  column  cross- 
section  a  range  of  loads  is  applied  to  establish  the  shape  and  the  peak  of  the 
column  load-deformation  curve.  This  curve  provides  the  data  for  the  column 
representation  methods  of  Section  2.1.9. 


The  strain  equation  for  the  temperature  application  step  for  the  column  is 
essentially  that  shown  by  Eq.  (2.1.16).  For  the  applied  load  step  in  the 
temperature-load  sequence  the  element  strains  are  found  by  the  following 
equation 


+  (ecip  +  Vi  +  (K°px+  Kwx+  Kp^j 


(2.1.48) 
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whtrt 


1  -  cos  *fl  /  e. 


m 


cos  ff/2  (eao/e( 


op  “cr- 


when  o  compression  fbocTIs  acting-. 

ep-  and  ^pxj  must  be  determined  by  an  Iterative  procedure  to  satisfy  the 
following  equilibrium  conditions. 


2  FM;  AM 
®ap.  =  M  *  .. 

1  *  emam 


^  EMAM 


(2.T.49) 


In  Eq.  (2.1.47)  the  stresses  F^.  are  determined  by  the  procedures  of  Section 
2.1.5.  The  iteration  equation:?  for  the  equilibrium  condition  are: 

,  ,  .  .  .£(FMi>r-lAM 

(epi>r=(epi)r-l  +  «opj-  - 

i  EMA« 


^pxj^r  ^pxj^r-1 1  ^apxj  +  ^^wxj^r-1 

_  fr>  \  A  . 


'  s.  (FM|)r-lAM  ^ 
fka-^ka  /^M\  ^ 


(2.1.50) 


fit 
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Two  sets  of  answers  are  provided  at  convergence  of  the  Iteration  equations 
(2.1.50)  with  the  second  set  defining  a  post-buckling  point  on  the  column 
foad-deformatfon-  carve. 

The  column  procedures  presented  In  this  section  are  based  on  constant 
moments  and  a  maximum  value  of  K^x]*  column  studies  have  shovel 

this  approximation. to  be  quite  accurate  for  prqctlcal  design  use.  Compar¬ 
isons.  of  column  food-deformation  curves  using  the  approximate  procedures 
of  this  section  and  by  Integrating  the  Inelastic  effects  a  Tong  the  length  of 
the  column  show  good  correlation  between  the  two  methods  with  the  approxi¬ 
mate  procedures  being  slightly  conservative.. 
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2.1.9  COLUMN  REPRESENTATION 


Since  many  structures  contain  column  members  in  the  form  of  stringers,  truss 
members,  etc.  it  is  most  convenient  to  consider  the  column  as  a  single 
element  in  the  cross-section  when  defining  the  allowable  lood-deformation 
curve  using  the  general  allowable  stress  procedures  described  ?n  Sections 
2.Y.?  through  2.1.7.  The  load-deformotior  characteristics  of  the  column 
may  not,  however,  be  defined  directly  by  a  simple  stress-strain  relationship 
such  as  the  Ramberg -Osgood  representation  of  material  property  stress-strain 
curves.  The  post-bockling  portion  of  the  column  curve  os  shown  i»r  Section- 
2.1.8  is  defined  by  bending  moment  equilibrium  where  the  applied  axial 
Toad  must  be  fess  than  the  critical  peak  buckling  load.  The  column,  however. 
Is  analogous  to  the  plate  buckling  problem  since  the  critical  buckling  strains 
of  both  are  based  on  the  Euler  equation.  Therefore,  a  good  approximation 
of  the  column  load-deformation  curve  can  be  obtained  by  the  use  of  the 
effective  area  coefficient  for  any  column  element  n  as  follows: 


(2.1.51) 


as  described  in  Sec.  2. 1 .6  of  this  report  and  Sec.  7-5  of  Ref.  (1).  This 
allows  on  equivalent  Ramherg-Osgood  stress-strain  curve  defined  by  Eq. 
(2.0.1)  to  be  used  in  the  calculations  for  the  particular  cross  section. 


The  equivalent  Ramberg-Osgood  stress-strain  curve  for  the  coiumn  element 
must  be  obtained  for  the  particular  temperature  distribution,  geometry,  and 
material  variations  in  the  column  cross-section  for  the  unrestrained  condition. 
From  this  stress-strain  curve  the  slope  representing  the  equivalent  modulus  of 
elasticity,  the  equivalent  yield  stress  defined  by  Esec  =  0.7E,  and  the 
Ramberg-Osgood  stress-strain  shape  parameter  m  are  obtained.  These  values 
are  used  in  the  structural  cross-section  to  define  the  equivalent  Ramberg- 
Osgood  stress-strain  curve  represented  by  Eq.  (2.0.1). 


For  the  column  element  in  the  cross-section  analysis  a  value  of  ecrn  must  be 
determined  to  define  the  effective  area  coefficient  given  by  Eq.  (2.1.32). 
First,  the  column  load-deformation  curve  is  computed  by  the  procedures  of 
Sec.  2.1.8  including  the  effects  of  local  buckling.  From  the  column  load- 
deformation  curves  the  two  significant  values  required  to  define  the  equiva¬ 
lent  critical  buckling  strain  are  the  peak  stress  defined  by  (Fap)  max,  and 
the  foreshortening  strain  at  which  the  peak  stress  occurs  defined  by  emax: 


The  equivalent  ecrn  for  the  column  element  is  defined  by  the  following 
relationships.  For  load  equilibrium  on  the  column  cross-section  at  the  strain 
associated  with  the  peak  stress: 
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2E«*M  P-rj2I 

where  the  subscript  M  denotes  element  number  In  the  column  cross-section, 
and  n  is  column  element  n  In  the  composite  structure.  In  terms  of  the 
effective  stress  on  the  total  cross-section  area  of  the  column 

Substituting  Eq.  (2.1.49)  IntoEq.  (2.1.51) 

(Fc)e  fecrn)*^  ^“^max.  SEMAM 

£AM 

Solving  for  ecrn 

(ecrn)1/2  =  Kp)ma^  (emax.)  *  (§)e 

XAM  (2.1.55) 


ecrn  (emax.) 

(CaP)max  SEMAm' 

2  =  (emax.) 

{FaP)max. 

(Fc)e  £AM 

• 

_«y. 

where  emax  is  the  foreshortening  strain  associated  with  the  peak  stress  of 
the  column  load-deformation  curve  and  (Fc)e  is  the  stress  associated  with 
emax.  as  determined  from  the  equivalent  Rambeig-Osgood  compression  stress- 
strain  curve. 

Since  each  element  in  the  column  cross-section  may  have  a  different  area, 
modulus  of  elasticity,  yield  stress,  temperature,  and  coefficient  of  thermal 
expansion,  an  equivalent  value  for  the  temperature  expansion  strain  is 
required  for  the  solution  of  more  complex  structures  which  treat  the  column 
as  a  single  structural  element.  This  temperature  expansion  strain  may  be 
expressed  as: 

(^Oequiv.  =  et  =  A^C^M  (2.1.56) 

2EMAM 
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where  the  subscript  M  denotes  element  number  Tn  the  cofunm  cross-section. 

Input  data  for  the  Allowable  Stress  IBM  Program  (E2I requires  separate 
input  values  for  a  and  T  for  each  element.  Therefore,  when  using  fhfr  prr?“ 
gram  for  cross-sections  which  contain  columns  as  a  single  element  it  may  be 
necessary  to  compute  an  equivalent  temperature  change  from  datum.  For  the 
column  as  a  single  element  in  a  cross-section  select  a  reference  volue  of  a  R 
and  compute  the  equivalent  temperature  change  by 

OJeqinv.  *  cMAmC«%  &T.ST 

°rSeMAM 

COcquiv.  aR  are  f^e  'nput  values  for  the  column  element. 

The  curves  on  page  40  demonstrate  the  use  of  the  effective  area  approxi¬ 
mation  and  provide  Justification  for  the  use  of  this  approximation. 
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2.1.10  FAILURE  CRITERIA 

Irr  crier  to  specify  thermal  limits  It  Is  necessary  ta  derive  suitable,  attawabl*. 
stresses  for  individual  elements  in  the  complex  structure.  Standard  allowables, 
such  as  for  crippling,  buckling,  and  compression  yield,  are  relatively  simple 
to  describe  for  variations  in  temperature.  Other  allowables,  such  as  honey¬ 
comb  general  buckling  are  more  difficult  to  acquire. 

In  determining  the  critical  honeycomb  panels  to  be  analyzed,  six  passible 
modes  of  failure  must  be  considered.  The  first  and  simplest  mode  is  a  tension 
failure  of  the  tension  face  due  to  bending  in  the  panel.  Another  failure, 
which  can  Be  produced  by  bending  is  core  shear  for  which  the  allowable  may 
be  expressed  asr 

qxcore  =  <*c/S)2  1/2  .2b  (2.1.58) 

The  loads  inducing  failure  in  this  mode  are  not  appreciably  changed  by  the 
thermal  loading  and  therefore  is  not  considered  as  on  important  failure  mode 
for  this  study.  A  third  failure  mode  in  bending  which  con  be  considered  is 
core  crushing,  the  allowable  for  which  is  expressed  ok 

Fc  =  (C]Fcy)  2/3  (KCE<)1/3  2(tQ/S)5/3  =  (2*F|)  FCF*> 

hEs  (2.1.59) 

This  mode  of  failure  has  been  found  not  to  occur  in  practical  design  cases 
and  is  therefore  also  not  considered  in  the  present  study.  A  fourth  mode  of 
failure  which  can  occur  in  either  bending  or  compression  panels  is  wrinkling 
of  the  facing  for  which  the  allowable  is: 

Fcw  =  0.5  jo.055  (tc/S)2  GCEC  (Es  +  3ETj  1/3  (2.1.40) 

The  fifth  failure  mode  which  may  occur  in  bending  os  well  as  compression  Is 
called  intercell-buckling  which  has  an  allowable  expressed  as; 

Fc.=0.9  Er  (tF]/S)  S/2  (2.1.61) 

The  honeycomb  allowable  Equations  2.1.60  and  2.1.61  can  be  computed 
independent  of  the  final  structural  arrangement  in  a  general  form  of  ecj  or 
Fcj  Vs  tp/S  and  eca  or  Fcw  VS  tc/S.  Such  a  plot  is  then  needed  only  for 
each  material  used.  Therefore,  to  reduce  the  computing  time  of  the 
"Allowable  Stress  Program"  this  data  is  entered  as  a  tcble  rather  than  being 
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computed  each  time  through  the  program.  Furthermore,  It  was  found  that  fay 
plotting  ec;  and  ecwpgalnst  the  respective  porameten  there  was  no  signi¬ 
ficant  change  with  temperature  within  the  realistic  range  of  the  allowable 
strains.  Thus,  if  the  calculation  of  data  is  done  for  a  medium  range  temp¬ 
erature  the  data  will  be  satisfactory  for  the  full  temperature  range  (see  Fig. 
2.3  and  2.4).  The  normal  allowable  strain  range  for  the  2024-T4  m ^terial 
plotted  In  Fig.  2.3  and  2.4  Is  from  .003  to  .0065.  Farther  discussion  of  the 
ISfA  70?  computer  program  for  producing  this  data  cao  be  found  In  eppendi« 
D  of  this  report. 


The  ffnaf  failure  mode  can  occur  only  fn  compression  of  the  panel  and  Is 
referred  to  as  general-buckling  for  which  the  allowable  U. expressed  me 


FccR  =  n2  *Fl  »F2  C  +  *Fl  +  *1? 

PM  \  -2  ) 


c2  0  -  U2)  Ofj  ES]  +  *P2  Es^  ^Fj  + 


(2.1.625 


i  +  L  +  Vx  Vy  F 
“04 


K  =  Cj  +  2C2  +  C3 

l  ■  (Ci  +  1-  c9)  vx  +  (o*+  \  -  n  c9)  vy 

2  C4  ~~i 

F  =  C|C3  -  C22  +  I  -  F  C2K 


Vx  =  2  C  tF]  tp2  E]E2 _ 

(»F|E|  +  tp2E2)  ^  G#cx  0  -M2) 

Vy  =  2  C  tFl  tp2  E]E2 _ _ 

(tp^j  +  tF2E2)  o2  Gey  (1  -  ft2) 

G>cx  =  effective  shear  modulus  of  core  in  x  direction 

G?cy  =  effective  shear  modulus  of  core  in  y  direction 

The  constants  C]/  C2>  C3/  C4  are  dependent  upon  the  edge  fixity  conditions 
of  the  panel.  A  listing  of  these  values  follows: 


42 


NORTH  AMERICAN  AVIATION  .  INC. 

COtVMlUS  »tvisto« 


Edges  simply  supported 

C*  *  *  G*  «  M 

73T-  “P- 

b  =  panel  length 

n  =  number  of  half  waves 

Loaded  edge  simply  supported 
Unloaded  edge  clamped 

Cl  =  5.33  b2  C2  =  1.33  C3  =  n2a2  C4=  1.33  b2 

b2 

Loaded  edge  clamped 
Unloaded  edge  simply  supported 

C2  =  1  C3  =7  n^  +  6n2  +  1»  a2 

‘V*'  ■’  ?■ 


For  n  =  1 

Ci = .75  b2 

T 

C4=3  b? 

7 

For  n  >1 

c,=c4=  1 

b2 

n2  + 

“  7 

AH  edges  clamped 

C2=  1.33 

C3  =  (n^  +  bn2  ■ 
n2  +  1 

f  1)  a2 

-  IT 

For  n  =  1 

Ci  =  4C,  =  4  b2 

.  1  4  -y 

0 

For  n  >1 

Ci  =4Q  =  5.33 

(n^ 

b2 

a 

These  last  three  failure  modes  are  the  most  important  in  the  design  of 
honeycomb  panels  and  constitute  the  primary  failure  modes  found  in  con¬ 
ventional  honeycomb  structures  subjected  to  thermal  stresses.  These  modes 
along  with  the  tension  mode  make  up  the  failure  modes  used  in  this  study. 
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Unlike  the  allowables  and  F,^,  Fccg  most  be  computed  within  the 
allowable  stress  program*  This  computation  takes  place  within  subroutine 
HKM  which  is  diagrammed  In  Appendix  D.l  of  this  report*  The.  Eqpatior» 
2*1.62  has  been  substantiated  by  a  test  program  presented  In  Appendix  A 
and  discussed  In  Section  4.0  of  this  report. 

Equations  (2.1.58)  through  (2.1.62)  were  substantiated  by  tests  performed 
by  the  NAA  Los  Angeles  Division  (See  Reference  3).  A  production  full 
depth,  honeycomb  component  was  tested  and  failed  In  the  made  of  face  . 
wrinkling,  showed  very  good  agreement  between  test  data  ond  the  appficaNfe 
allowable  stress  Equation  (2.1.60)  (Ref.  4,  Figure  D.5). 

For  tension  elements  in  the  structure  the  failure  criteria  is  defined  as.  a 
maximum  strain  or  elongation  of  the  material.  Compression  elements  which 
buckle  locally  will  reach  a  peak  load  causing  further  increases  of  load  and/ 
or  moment  to  be  taken  by  remaining  members  of  the  structure  when  some 
member  will  eventually  reach  the  strain  cut-off  or  the  crippling  cut-off 
defined  by  the  following  equation. 

Crippling  of  compression  elements  is  considered  to  occur  if: 

FnCn  >  Fyn  (2.1.1 

A  tension  ultimate  stress  cut-off  Is  used  on  all  Ramberg-Osgood  stress-strain 
curves  with  a  compression  yield  stress  cut-off  used  on  all  Ramberg-Osgood 
stress-strain  curves  used  to  represent  a  column  element. 
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2.1.11  LOAD-DEFORMATION  CURVES 

The  previous,  wctims  describe  the  procedures  for  obtaining.  the  «l>— t 
strains  and  stresses  at  any  particular  step  In  the  temperature^Toad  sequence. 
Convergence  of  the  iteration  procedures  of  Section  2.1.7  produces  the  final 
values  of  ep •,  KpX.,  and  K  .  at  each  step  in  the  temperature-load  sequence. 

At  any  step  \,  the  tatal  inelaAlc  effects  for  all  steps  on  any  element  n  in  the 
cross-section  may  be  expressed  by 

ePsn  =2j  ep.  +  (£j  yn  +  CS.  Kp^)  x„  (2.1.6^ 

where  ep$n  represents  the  permanent  set  or*  each  element  of  the  cross  .section 
if  step  |  ends  with  the  removal  of  all  load  ctnd  temperature  and  assuming  that 
the  residual  stresses  do  not  change  due  to  creep  at  room  temperature. 

If  the  procedures  of  the  previous  sections  are  followed  forgiven  geometry, 
materials,  temperature  distribution,  and  temperature-load  sequence  for  a 
range  of  values  of  applied  axial  loads  and  bending  moments  then  an  allow¬ 
able  load  or  moment  vs.  strain  curve  can  be  constructed.  In  this  case,  the 

primary  bending  moment  (Mj^  may  be  plotted. against  Fapn^Em  +  ep$  •  If 

is  essential,  however,  that  fixed  ratios  exist  between  Pj,  M*.,  ana  Myj  if 
a  single  moment-strain  curve  is  to  be  used  to  define  the  first  instability 
failure,  yield  and  ultimate  load  for  the  cross-section. 

For  a  range  of  values  of  load  and/or  bending  moment,  the  elastic  strain 
term  (Fapn/f^  may  be  defined  as  follows: 

FaP«/Em  =  *1  e°pi  +  (2i  KoPx)Z2L  +  (si  ’S'j)  {2  U65) 

where  m  designates  the  element  with  the  maximum  permanent  set  as  defined 
by  Eq.  (2.1.64). 
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2.2  SAMPLE  PROBtEM 

To  provide  substantiation  for  the  newly  developed  methods  presented  In  this 
report  a  sample  problem  on  the  computer  was  made  for  comparison  with  tests 
performed  during  this  study.  Two  comparisons  are  made  in  this  section.  The 
first  comparison  and  computer  run  was  made  to  predict  the  allowable  load  for 
the  skin  stringer  columns  found  in  the  test  horizontal  stabilizer.  'Hie  corn-  » 
pu ter  run  was  made  with  the  column  allowable  program  E2140  described  br^ 
Appendix  C.  1 .  This  column  allowable  was  computed  without  the  lateral 
bending  effect  discussed  in  Section  2.2.1.  The  lateral  bending  was  not  used, 
in  this  case  because  of  the  test  conditions  would  not  produce  the  tension 
loads  caused  by  buckling  because  of  the  lack  of  edge  restraint  (see  Figure 
B.30).  Also  because  of  the  end  restraint  conditions  the  column  end  fixity 
was  a  different  value  in  this  comparison  problem  than  in  the  test  horizontal 
stabilizer  problem.  The  results  of  this  computed  column  allowable  and  the  ■ 
test  allowable  can  be  found  in  Section  2.2.3. 

The  second  comparison  and  computer  run  was  made  to  predict  the  ultimate 
failure  load  for  the  typical  horizohtal  stabilizer  with  an  induced  temperature 
gradient..  Before  making  the  actual  computer  run  it  was  necessary  to  compute 
several  groups  of  input  data.  The  first  group  of  data  required  was  the  honey¬ 
comb  allowable  strains  against  the  parameters  tp/S  and  tg/S  as  a  function  of 
temperature.  This  data  was  compiled  through  the  use  of  the  computer  pro¬ 
gram  E2128  described  in  Appendix  C.2  of  this  report.  This  data  is  then  used 
as  a  table  input  as  described  in  Appendix  D.1.3  for  the  final  computer  run. 

The  second  group  of  data  required  is  compiled  through  the  use  of  computer 
program  E2140  described  in  Appendix  C.  I.  In  addition  to  the  straight  axial 
load  a  laterally  induced  bending  moment  was  also  used  as  discussed  in 
Section  2.2.1.  This  data  was  then  used  as  element  buckling  coefficients 
for  the  Input  data  such  as  described  in  Appendix  D.1.3. 

Another  important  part  of  the  input  data  is  the  material  properties  E,  Fty, 
Ftn,  and  the  Ramberg-Osgood  coefficient  (m).  The  derivation  of  this  data 
for  the  specific  test  horizontal  stabilizer  problem  is  discussed  in  detail  in 
Section  2.2.2.  Although  tensile  coupon  tests  were  made  in  this  case  to 
Insure  greatest  accuracy  if  was  found  that  the  use  of  Reference  5  would  give 
very  satisfactory  results. 
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2.2.  T  SKI N-STRl MGER  COMPRESSION  ALLOWABLES  INCLUDING 

THE  EFFECTS  OF  INDUCED  LATERAL  BENDING 

The  allowable  compression  stresses  (Fcr)  or  the  critical  buckling  strnihr 
(eCr)  of  skin-stringer  column  elements  may  not  be  accurately  determined  by 
the  application  of  direct  compression  loads  either  analytically  or  experi¬ 
mentally  since  the  compression  load  alone  does  not  always  represent  the 
only  load  acting  on  the  stringer.  More  representative  allowables  could  be 
determined  experimentally,  however,  by  box  beams  with  applied  bending- 
moments  where  the  chordwise  restraint  to  the  skin  panels  and  the  effectr 
of  radius  of  curvature  are  present.  Essentially,  the  major  effect  of  the  bend¬ 
ing  radius  of  curvature  of  the  Box  beam  is  to  force  inward  buckling  of  the 
skin  panels  between  the  stringers.  This  inward  buckling  on  both  sides  of  a 
stringer  coupled  with  chordwise  restraint  of  the  skin  panels  produces  chord- 
wise  tension  stresses.  Components  of  these  stresses  normal  to  the  plone  of 
the  stringer  produce  bending  moments  on  the  stringer  which  are  acting  simul¬ 
taneously  with  the  applied  compression  loads.  In  many  structures,  these 
.  effects  may  be  quite  small  but  the  induced  bending  moments  In  shallow  depth 
stringers  may  produce  substantial  losses  in  allowable  compression  strength. 

The  skin-stringer  panels  in  the  typical  test  horizontal  stabilizer  represent 
the  type  of  structure  where  the  induced  lateral  bending  due  to  skin  panel 
buckling  should  be  investigated.  The  following  analysis  Is  a  conservative 
approximation  to  allow  for  the  induced  bending  effects. 

One  typical  skin-stringer  pane!  is  considered  where  bj  =  3.73  inches 
(between  stringer  £  end  ts  =  .042  in.  Temperature  is  uniform  at  375°E. 

For  the  panel  simply  supported  on  all  four  edges,  the  critical  buckling  strain 

ecr  =  K^t_J2  =  3-62  ^  =  .000459  in./tn. 

Considering  the  center  strip  of  the  buckled  plate  to  act  as  a  column,  the 
deflection  may  be  approximated  by 

Wm  =  2  P 

where  P  is  the  radius  of  gyration  of  the  strip  element,  ecr  is  the  buckling 
strain  of  the  panel,  and  (  A  L/L)  is  the  foreshortening  of  the  panel  after 
buckling.  This  foreshortening  is  taken  as  the  difference  between  the  strain 
at  which  the  skin-stringer  column  peak  occurs  and  the  critical  buckling 
strain.  From  Figure  2.5,  the  column  peak  occurs  at  .00365  in./in. 
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The  maximum  bending  moment  at  the  center  of  the  column  span  is 

M  =  wl?  =  59.7(8.17)?  =  497.5 In.  -tbs. 

8  8 

An  applied  bending  moment  rotational  strain  on  the  column  cross-section  is 
calculated  by 

-  Mt-  =  497.5  (.324  “•  -JB22&  Win,  at  mOmm » 
fiber.  “(1,)"“  73700 

The  values  of  e,^  for  the  column  elements  are  determined  by  the  procedures 
of  Section  2.1.9  with  some  typical  results  shown  in  Figti.e  2.5  for  the  critical 
column  element.  From  the  column  load-deformation  curves  ecrn  **  calcu¬ 
lated  in  the  table  below  for  each  stringer  in  the  compression  cover. 


COMPRESSION  STRESS,  Fc,  psl.  x  10*3 


STRAIN  ,  IN./IN 


TABLE  2.2 

Effective  Column  ecr  Value*  Including  the  Effect*  of  Induced  Latercj|  Bentjlno 
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2*2.2  MATERIAL  PROPERTIES' 


The  tension  stress-stroin  curves  used  In  the  analysis  and  represented  by  the 
Ramberg-Osgood  equation 


eE 


=  F 


are  shov-n  in  Figure  2.6.  The  curve  of  Figure  2.7  shows  the  data  from  two 
tensile  coupons  which  was  used  to  determine  the  Rambetg-Qsgood  shape, 
parameter,  m.  The  stress-strain  curves  shown  in  Figure  2.7  were  obtained 
from  standard  NAA  tensile  coupons  taken  from  the  test  stabilizer  between 
H.S.  Sta.  167  and  H.S.  Sta.  177. 


The  material  properties  obtained  from  these  coupons  showed  a  distinct 
permanent  loss  due  to  the  elevated  temperature  curing  of  the  high  temperature 
strain  gages  in  the  test  area  (H.S.  Sta.  187  to  H.S.  Sta.  197).  The  entire 
horizontal  stabilizer  was  placed  in  an  oven  and  heated  to  350°F  for  two 
hours  and  then  to  450°F  for  one-half  hour  for  the  curing  process.  The 
tensile  coupons  were  taken  from  the  stabilizer  skins  after  the  failing  ioad 
test  was  completed.  Therefore,  the  coupons  and  the  test  area  had  the  same 
temperature  history  prior  to  the  beginning  of  the  tests.  The  primary  cause 
of  the  permanent  loss  of  material  properties  was  the  temperature  of  450°F 
for  one-half  hour.  Oata  obtansd  from  short  time  elevdted  temperature 
tensile  coupon  tests  is  plotted  in  Figure  2.8  and  used  as  material  property 
input  data  for  the  analytical  portion  of  the  program.  The  room  temperature 
values  of  Fty  and  Fty  on  Figure  2.9  are  specifically  noted  for  comparison 
wi  th  the  expected  recovery  values  after  exposure  to  the  elevated  temper¬ 
ature  environment  required  for  strain  gage  curing.  Using  data  from  Refer¬ 
ence  (5)  and  Figure  3.2.7. 1.1  the  expected  room  temperature  recovery 
points  are  also  plotted.  From  Figure  2.9  the  room  temperature  recovery 
factors  are  found  to  be 

/FfuR.T.  =  °-785  ”d  FV/Vt.  =  °-m 


The  room  temperature  materia!  allowables  are  taken  as  Fty  =  78000  psi. 
and  Fty  =  68000  ps  .  from  Table  3. 2. 7.0  (b)  of  Reference  (5).  For  the 
temperature  time  history  of  the  stabilizer  the  room  temperature  material 
allowables  may  be  expected  to  recover  to  Ftu  =  0.785  (78000)  =  61200  psi. 
and  Ft  =  0.722  (68000)  =  49100  psi.  Although  good  correlation  is  ob- 
tained'from  room  temperature  recovery  properties  it  is  impossible  to  predict 
what  the  expected  recovery  values  should  be  at  temperatures  between  room 
temperature  and  the  maximum  soaking  temperature  due  to  a  lack  of  test 
data.  Since  it  appears  that  the  recovery  problem  can  be  quite  serious  it 
would  also  indicate  that  a  worthwhile  program  could  be  established  to 
determine  re~overy  properties  ot  temperatures  above  room  temperature  for 
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FIGURE  2.6 

tension  stress-strain  curves 

RECOVERY  CURVES  AFTER  EXPOSURE  TO  350°  F  FOR  2  HR 
AND  450°  F  FOR  1/2  HR. 
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FIGURE  2.7 

ROOM  TEMPERATURE  TENSILE  COUPON  STRESS-STRAIN 
CURVES  AND  EQUIVALENT  RAMBERG-OSGOOD  EQUATION 


MODULUS  OF  ELASTICITY,  E,  p* f.  x  10“ 
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EFFECT  OF  TEMPERATURE  ON  THE  MECHANICAL  PROPERTIES 
OF  7075-T6  BARE  ALUMINUM  ALLOY  FOR  A  TEMPERATURE  HISTORY 
OF  350°F  FOR  TWO  HOURS  AND  450°F  FOR  1/2  HOUR 
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*  Fty  IS  THE  0.002  IN./IN.  OFFSET  YIELD  STRESS 
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FIGURE  2.9 


MATERIAL  RECOVERY  PROfSTIK  . 


S.T.  =  78000  ^ 


SOAKING  TIME  AT  TEMPERATURE,  HRS 


PREVIOUS  TEMPERATURE  HISTORY 
2  HOURS  AT  350°F 
4  HOURS  AT  450°F 
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2-2  .3  COMPUTED  AND  TEST  DATA  COMPARISON 

A  test  was  performed  on  two  sfcm  stringer  panels  taken  from  the  test  horf- 
zontai  stabilizer.  This  test  was  made  to  substantiate  the  prediction  of  columt 
allowables  to  be  used  In  the  allowable  stress  computer  program.  The  test 
set  up  and  results  are  presented  in  Appendix  B.  The  test  specimen  was  failed 
under  a  varying  temperature  which  produced  the  desired  temperature  gradients 
Since  the  computer  program  cannot  exactly  simulate  the  varying  temperature^ 
the  allowable  food  was  computed  for  the  temperature  disfritution.in  the 
panel  at  the  time  of  failure.  The  comparison  of  test  and  computed  results  are 
shown  in  Figure  2.10.  This  comparison  shows  excellent  agreement  between 
the  maximum  test  load  and  the  load  computed  with  the  temperature  distribu¬ 
tion  at  the  time  of  failure.  From  these  tests,  it  was  concluded  that  the 
column  allowable  program  Appendix  C  was  satisfactory  for  estimating  column 
input  data  for  the  main  program. 

The  test  set  up  and  results  of  the  test  horizontal  stabilizer  are  presented  in 
Appendix  B.  Again  the  temperature  distribution  at  the  time  of  failure  was 
used  to  compute  the  allowable  moment  for  the  test  section.  The  test  section 
was  broken  into  37  structural  elements  os  shown  in  Figure  2.11.  The  honey¬ 
comb  element  allowables  (Fcw  Fc;)  were  computed  using  the  IBM  program 
E2128-presented  in  Appendix  cf.2.  The  skin  stringer  element  buckling  co¬ 
efficients  required  were  computed  from  the  column  allowable  program  E2140 
in  Appendix  C.l.  Allowables  for  other  elements  were  computed  as  described 
In  Section  2.1. 

The  first  computer  run  results  plotted  in  Figure  2. 12  showed  an  allowable 
bending  moment  Mx  of  186,000  in*.  This  moment  is  112%  of  the  actual 
test  failure  load.  At  this  point,  further  investigation  into  laterally  induced 
bending  loads  on  the  skin  stringer  elements  was  done  as  presented  in 
Section  2.1.8.  Following  this  investigation,  the  second  computer  nm 
plotted  in  Figure  2.12  was  made  showing  an  ultimate  allowable  bending 
moment  of  the  section  os  171,000  in*.  This  computed  allowable  moment  is 
103%  of  actual  failing  load. 

Figure  2.13  shows  a  comparison  of  measured  and  calculated  strains  vs.  chord 
at  H.S.  Sta.  192.  The  measured  strains  are  shown  for  an  applied  load  of 
4,390  lbs.  and  the  calculated  strains  ore  shown  for  the  nearest  calculated 
applied  load  of  4,410  lbs.  The  calculated  strains  are  shown  for  the  temp- 
eroture-foad  sequence  +T,  +P  and  reflect  the  effects  of  the  non-uniform 
temperature  distribution  as  well  os  the  applied  load. 

From  this  data,  it  has  been  concluded  that  the  proposed  allowable  stress 
method  developed  in  this  study  will  meet  the  requirements  for  predicting 
nuclear  effects  temperature  limits  within  the  accuracies  required. 


FIGURE  2.11 
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FIGURE  2. 12 

MOMENTOffORMATlOM  CURVE 
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CHORD  LENGTH 


COMPRESSION  SURFACE  AT  350°F  FOR  ZERO  LOAD 
TENSION  SURFACE  FOR  ZERO  LOAD 
COMPRESSION  SURFACE  AT  350°F  FOR  140%  (4390#) 
TENSION  SURFACE  FOR  140%  (4390#) 
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3J)  A  SIMPLIFIED  HAND  CALCULATION  METHOD 

The  hand  calculation  method  presented  In  this  report  was  selected  to  provide 
reasonable  accuracy  without  excessive  detail.  This  detail  which  accounts 
for  the  last  increments  of  accuracy  can  be  eliminated  from  the  thermal  stress 
analysis  to  greatly  simplify  the  computations  required. 

The  first  step  in  the  elimination  of  detail  was  made  by  eliminating  chord- 
wise  bending  effects  from  the  analysis.  This  chordwise  effect  accounts  for  a 
maximum  of  10%  of  the  total  stress  in  most  cases.  The  location  of  the  critical 
area  forward  and  oft  on  the  surface  actually  determines  the  magnitude  of  the 
induced  error.  If  the  critical  area  being  analyzed  is  close  to  the  forward 
and  aft  center  of  area,  essentially  no  errors  exist.  If  the  critical  areas 
ere  on  the  leading  or  trailing  edge  of  the  surface,  an  error  of  +  10%  may  exist. 
The  existence  of  this  error  is  dependent  upon  the  temperature  distribution  on 
the  structure.  An  evaluation  of  the  thermally  induced  drag  moment  will 
indicate  the  sign  of  the  error. 

A  direct  elimination  of  detail  can  be  made  in  the  selection  of  structural 
elements  to  be  included  in  the  analysis.  The  object  of  element  selection  is 
to  group  similar  detailed  elements  according  to  failure  mode,  temperature, 
and  distance  from  the  x  axis.  Using  this  procedure,  all  similar  column  ele- 
mer  ts  with  approximately  the  same  temperature  and  moment  arm  can  be  con¬ 
sidered  as  one  element  at  an  average  temperature,  moment  arm  and  buckling 
allowable.  The  induced  error  from  such  a  procedure  should  remain  within 
+  2%.  The  magnitude  of  this  error  is  dependent  upon  the  accuracies  with  which 
the  average  values  are  computed  for  the  combined  elements.  Selection  of 
typical  elements  Is  demonstrated  in  the  sample  problem  of  this  section. 
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3.T  SIMPUFIEP  METHOD  DEVELOPMENT 
3.1.1  THERMAL  STRESS  EQUATION 


Thermal  stress  distributions  can  be  calculated  for  appropriate  flight  loads  and 
temperature  distributions  with  a  general  equation  (3.1.7)  which  was  derived 
by  an  expansion  of  information  obtained  from- Reference  (I) as  folio* 


Assume  a  one-dimensional  temperature  distribution  T(y). 


With  the  assumed  temperature  distribution  T(y)  and  the  plate  fully  restrained 
with  no  buckling,  the  stress  distribution  becomes: 


fx  *  “ETW 

«  =  coefficient  of  thermal  expansion  (3.1.1) 

E  =  modulus  of  elasticity 


If  we  free  the  plate  of  end  but  not  bending  restraints  a  relief  stress  is  realized. 
This  relief  stress  is  superposed  on  the  full  restraint  stresses.  The  total  tensile 
load  Px  may  be  written  as: 


«ET(y)  bdy 


The  resultant  relief  stress  becomes 

(gt=_Px=j^L  f  «ET(y)dAt  (3.1.2) 

Aj  At 

A|  =  total  area  of  section 
dAt  =  b  dy 


Since  the  temperature  distribution  is  unsyrnmetricol  the  tensile  forces  have 
c  resultant  moment.  Thus  when  the  plate  is  unrestrained  in  bending  an 
unbalanced  relief  moment  results: 


66 


NORTH  AMERICAN  AVIATION  ,  INC. 


COtUMtut  OrviVON 


Mr'  r 


£T(y)  bydjr 


The  resulting  relief  stress  becomes: 


(ULu  =  -  MzV=.  Y  f  *EJ(y)ydAt 

vr- 


C*.U3) 


Y  =  distance,  to  element  C.G- 


The  one  dimensional  thermal  stress  equation  is  formed  by  compiling  Equations 
3.1.1,  3.1.2,  and  3.1.3  os  follows: 

f  =  aET(y)  ■  1_  f  a  ET(y>  dA.  -  J_  r  aET(y)ydA,  (3.1.4) 

Aj  '*  J. 


The  stress  from  the  external  bending  moment  Mz  can  be  expressed  as: 


fx  =  MY 


(3.1.5) 


Thus  by  combining  equations  3.1.4  and  3.1.5,  the  final  stress  equation 
becomes: 

fx  = -aET(y,z)  +  1  r  aET(y)dA^+  y  C  aET(y)dAj  . 

*^t  i  !z  *c 


+  Mzy 


(3.1.6) 


For  the  use  of  numerical  integration  Equation  3.1.6  becomes: 


fx  =  £ ET(y ,  z)_  Jl V  a  EAT(y)-  y  V  aEAyT(y)  +  J^y  (3.1.7) 

a  .  i 


where: 


A  =  area  of  small  elements  in  section 


Equation  3.1.7  is  modified  for  stresses  above  the  proportional  limit  by  the  use 
of  the  Ramberg-Osgood  Equation  as  follows: 
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E  becomes  Eg  *  . —  0,1 

1  +  3/7  {F/F^^ 

F  =  actual  stress 
Fy  =  yield  stress 

m  =  coefficient  dependent  on  materia!  stress  strain  curve 

2  equation  fs  also  modified  for  buckling  by  changing  the  effective  area 
a  buckling  element.  The  follovring  equation  is  used  to  modify  Ai 

AEFF  =  A  /Fc|/P  (3.1 

FcR  =  buckling  stress  of  element 
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3.1.2  ANALYSIS  DEVELOPMENT 

The  development  of  the  thermaf  stress  equation  Into  a  usoble  thermal  limits 
analysis  requires  several  complementary  steps.  Each  of  these  steps  are  pre¬ 
sented  here  In  a  form  to  simplify  the  computations  as  much  as  possible.  The 
developed  analysts  method  ts  also  readily-  adaptable  fa  small  digital  qc 
type  calculators. 

The  first  step  tn  deriving  a  thermaf  limit  for  a  complex  aircraft  structure  ft- 
to  determine  the  temperature  distribution  on  the  structure.  Knowing  the  Mach 
numi  ers  and  altitude  at  which  the  limit  is  desired  a  curve  of  skin  thickness 
vs  average  element  temperature  for  the  given  boundary  layer  temperature 
and  reference  AT  rise  con  be  constructed  such  js  Figure  3.2.  From  such  a 
curve  an  approximate  temperature  distribution  can  be  determined  for  any 
given  structure.  It  is  recommended  that  a  reference  AT  rise  be  selected  in 
the  range  of  from  300  to  400°F  in  .032  aluminum  to  insure  some  proximity 
to  the  final  limit  rise.  The  respective  honeycomb  element  temperatures  must 
be  computed  for  the  specific  configuration  to  complete  the  structural  temp¬ 
erature  distribution. 

At  this  time  the  structure  should  be  broken  down  into  structural  elements  as 
described  in  paragraph  3.0  of  this  report.  With  the  number  of  elements 
end  temperatures  determined  a  table  of  computations  con  be  started  to  de¬ 
termine  the  stress  distribution.  The  following  is  a  listing  by  column  of  the 
necessary  tabular  data  and  computations  to  arrive  at  an  approximate  Stress 
distribution. 
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Column  1.  Element  number  or  ideoii RcsriMft 

2.  The  actual  area  of  each  element  of  Column  1. 

3.  The  temperature  rise  of  the  element  above  room  temperature. 

4.  Young's  Modulus  Ee  at  the  element  temperature. 

5.  Corrected  Area  computed  by  the  following  equation; 


A nax%/V& 


(3.1.W 


The  moment  arm  Y  of  each  element  from  a  reference.  X  axis.. 
The  products  of  the  Corrected  Area  and  moment  arm  Y  for 
each  element  are  summed  and  divided  by  the  summation  of 
column  5  to  compute  Y. 

The  corrected  moment  arm  Ycorr  listed  for  each  element, 
computed  by  the  following  equation; 


Yconr  V  -  Y 


(3.1.1 


The  moment  of  inertia  of  each  element  about  the  X  axis 
summed  to  give  the  totol  moment  of  inertia  of  the  section. 
Computed  by  the  following  equation: 


Ixx  -  Aeon-  x  Y 


(3.1.1 


The  full  fixity  thermal  stress  computed  by  equation  3. .1.1. 
The  summation  of  the  thermal  tension  relief  load  (£5  x  10), 
and  the  thermal  relief  stress  as  computed  by  equation  3.1.2 
where  Aj.  =  Acorr. 

The  summation  of  the  thermal  bending  relief  moment  (8  x  ii) 
as  computed  by  equation  3.1.3. 

The  applied  bending  stress  as  computed  by  equation  3.1.5 
(the  bending  moment  used  should  be  for  the  normal  load 
factor  Nz  of  greatest  interest). 

Iz  =  Ixx  o*1'*  Y  =  YCorr* 

Thermal  bending  relief  stress  computed  by  using  the  equation 


*b»h  "  ^th  -Ycorr 


(3.1.1 


Where  Mfh  -  ^Column  12. 
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Column  15.  Total  thermal  stress  on  each  element  as  computed  by 
Equation  3.7.4. 

16.  The  actual  corrected  applied  bending  stress  on  each  element 
os  computed  by  the  following  equation: 


fbact  “  *T>  \orr 

A 


(3.1.14) 


17.  The  actual  corrected  total  thermal  stress  on  each  element 
os.  computed  by  the  following  equations 


*th  oct  ~  *th  ^corr 
A 

18.  Final  total  element  stress  (16+  17) 


(3.1.15) 


Upon  the  completion  of  the  computation  of  the  Final  Total  element  stress 
each  element  stress  must  be  checked  against  the  respective  element  buck¬ 
ling  allowable  or  yield  stress,  whichever  is  applicable.,  with  the  buckling, 
if  critical,  taking  precedence.  If  an  element  Is  found  to  be  stressed  above 
the  allowable  buckling  a  correction  must  be  made  to  the  effective  area, 
Acorr/  using  equation  3.1.9  where  A  is  the  last  Acorr  and  Appp  is  the  new 
Acorr-  If  or>  element  stress  is  found  which  exceeds  .6  Fey  a  correction  of 
the  Acorr  must  be  made  for  this  element  also  using  the  following  equation: 


Area  =  (NEW)Acorr 

(3.1.16) 


Where  Es  is  computed  using  equation  3.1.8  with  E  =  EQ  and  F  is  the  final 
total  element  stress.  Upon  correcting  the  areas  for  these  elements  the  com¬ 
putation  must  be  again  performed  and  repeated  until  the  final  total  element 
stresses  converge  within  required  error  limits. 


When  final  element  stress  convergence  has  been  accomplished  the  element 
nearest  to  becoming  critical  or  most  critical  can  be  selected.  The  approx¬ 
imate  limiting  temperature  rise  or  moment  can  then  be  estimated  by  one  of 
the  following  methods: 
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3.2  SAMPLE  PROBLEM' 

The  sample  problem  computed  in  tables  3.)  through  3.5  were  computed  far 
a  typical  horizontal  stabilizer  at  a  flight  condition  of  Mach  1.2  at  20,000 
ft.  altitude  and  a  reference  temperature  rise  of  450°F.  on  .032  aluminum 
skin.  Figure  3.3  shows  a  layout  of  the  horizontal  stabilizer  section  and 
areas  included  in  each  structural  element  used  in  this  analysis. 

The  most  critical  element  in  this  analysis  is  element  9  which  is.  the.  full 
depth  honeycomb  trailing  edge.  The  final  element  stress  was  computed  to 
be  31/800  psi.  The  allov/able  stress  for  this  element  at  483°F.  is  computed 
using  Equation  (3.1.2),  the  reference  rise  orr  -.032 skin  becomes: 

450°  =  reference  rise  used  in  table 

450°  (  33500  ~  4100  )  =  476°F  allowable  rise  on  .032  at  N7  =  2 
V  27700  ' 

If  the  change  in  temperature  rise  was  great  enough  in  this  problem  a  correction 
of  the  allowable  for  the  increased  or  decreased  temperature  should  be  accounted 
for  by  recalculating  the  allowable.  The  allowable  rise  computed  here  agrees 
reasonably  well  with  the  allowable  calculated  for  the  same  condition  using  past 
methods.  The  single  mission  limit  for  this  condition  using  the  more  exact 
method  was  computed  as  44 5°F  rise  on  .032  during  that  study.  This  shows  a 
6.3%  error  in  the  proposed  hand  calculation  method  for  this  problem. 
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If  the  element  thermo!  stress  (column  16)  Is  greater  than  the  element  allow¬ 
able  stress  the  foi  Jawing  equation  should  be  used: 


(Tele  -  TBOUNDARYXFaHow  -  fb  act)  =  (TELE  -  TbOUNDARY) 
.  fth  DC- 


allow 

(3.1.17) 


fj^act  =  column  17 

Fjf,  oet*=  column  )6 

^ailow  =  element  allowable  stress 

TELE  =  element  total  temperature 

^BOUNDARY  =  Boundary  Layer  temperature 

The  new  allowable  element  temperature  rise  calculated  is  for  the  moment 
used  in  the  stress  distribution  table.  From  this  moment  a  value  of  Nz  can  be 
obtained  to  plot  versus  the  limit  temperature  rise.  By  changing  the  applied 
moment  and  fj,  in  the  table  proportionately  another  approximate  limit  temp¬ 
erature  rise  can  be  obtained  without  recomputing  the  entire  table. 

CASE  il 

If  the  thermal  stress  is  less  than  the  allowable  and  '/he  total  stress  is  greater 
than  the  allowable  the  applied  bending  moment  can  be  reduced  by  the 
following  equation: 

M(^ailow  -  fyh  act)  =  fallow  -(3.1.18) 

act 

This  allowable  moment  is  for  the  temperature  rise  used  in  the  computation 
of  the  actual  stress  distribution.  If  the  allowable  moment  should  become  less 
than  that  at  an  Nz  =  1  the  limit  temperature  rise  should  be  adjusted  as  in 
Cose  I. 

CASE  II! 

If  the  allowable  stress  is  greater  than  the  total  stress  either  the  moment  or  the 
temperature  rise  limit  can  be  computed  using  the  respective  equations  3.1.17 
end  3.1.13  cf  the  preceding  paragraphs.  1°  cny  of  the  discussed  equations 
if  the  element  temperature  is  changed  the  ollowoble  element  stress  should 
be  cdjusted  accordingly  thus  requiring  sc ~<r  iteration  in  the  computation  of 
an  allowable  temperature  rise. 
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4.0  GENERAL  HONEYCOMB  PANEL  BUCKUNG 

The  general  bockITng  equation  was  derived  Tn  Reference  6.  The  equation 
was  derived  to  handle  sandwich  panels  with  faces  of  unequal  thickness  and 
dissimilar  materials.  This  equation  then  fits  the  general  production  honey¬ 
comb  panel  exposed  to  thermal  radiation  exactly.  The  production  honey¬ 
comb  panel,  with  faces  of  unequal  thickness  because  of  external  damage 
requirements,  have  been  difficult  to  handle.  The  addition  of  a  temper  close 
gradient  and  change  of  material  properties  especially  in  the  hot  facing  has 
complicated  the  problem  to  a  point  where  testing  of  panels  was  the  only 
sofotioiT.  While  no w,  through  the  use  of  the  generalized  panel  buckling 
equation,  the  differences  of  facing  thickness  and  material  properties  can  be 
accounted  for  theoretically.  With  this  capability,  aircraft  thermal  limits 
can  be  derived  directly  through  the  use  of  a  theoretical  analysis.  use 
of  the  general  buckling  equation  also  provides  the  capability  of  designing 
honeycomb  structures  to  withstand  the  anticipated  thermal  loading. 

The  testing  of  the  honeycomb  panels  and  configurations  are  described  in 
appendix  A.  The  panels  were  tested  to  provide  correlation  data  for  the 
general  honeycomb  buckling  equation  proposed  for  use  in  this  study.  The 
data  compiled  in  appendix  A  showed  very  good  repeatability  and  consistency 
from  one  test  to  the  other.  In  a  comparison  of  the  test  data  with  computed 
allowables,  agreement  was  wide  spread.  Portions  of  this  spread  con  be 
explained  in  several  ways.  The  allowable  equation  shows  a  high  sensitivity 
of  the  allowable  load  to  the  hot  face  temperature.  A  study  of  the  equation 
sensitivity  to  temperature  shows  that  this  sensitivity  increases  with  temp¬ 
erature.  See  Fig.  (4.1).  Reading  and  measurement  errors  for  the  thermo¬ 
couple  installation  used  on  these  tests  was  estimated  at/25  °F  and  -5°F. 

A  second  source  of  error  in  the  computed  allowable  results  from  the  fact  that 
a  temperature  gradient  exists  across  the  panel  to  an  extent  of  5%  lower 
temperature  rise  at  the  edge  as  compared  to  that  at  the  center  of  the  panel . 
This  results  in  a  lower  average  temperature  than  that  measured. 

An  evaluation  of  panel  construction  tolerances  was  also  made  to  establish 
induced  computed  allowable  errors.  This  study  revealed  that  an  increase 
of  the  core  depth  of  /  5%  had  a  direct  effect  on  the  allowable  by  / 5%. 

A  /  5%  tolerance  on  the  thickness  of  the  hot  facfe  produced  a  /  4%  change 
inlhe  panel  allowable  load.  A  thickness  change  in  the  cold  facing  produced 
negectable  allowable  load  change. 

The  test  data  and  computed  allowable  data  is  compiled  in  table  (4.1)  and 
plotted  in  fig.  (4.2).  The  computed  allowables  were  computed  using 
equation  (2.1.62)  which  was  also  used  in  the  “Allowable  Stress  Program" 
as  subroutine  HKN.  Although  the  comparison  of  the  computed  data  and  test 
data  shows  a  wide  range  of  scatter,  the  inaccuracies  of  the  test  and  the 
sensitivity  of  the  equation  must  be  considered  in  the  evaluation  of  the  equa¬ 
tion.  Since  the  equation  shows  good  results  compared  to  test  data  at  room 
temperature,  and  again  crosses  the  0%  error  line  in  the  area  of  600°F.,  it  is 
the  opinion  of  the  writer  that  the  equotlon  is  satisfactory  for  predicting 
thermal  limits  in  honeycomb  structures. 
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5.0  THERMAL  STRESS  SUBSTANTIATION  TEST 

A  test  program  was  conduced  in  the  Structures  Laboratory,  Columbus 
Division  of  North  Amencon  Aviation,  Inc.  in  support  of  a  Generalized 
Thermal  Stress  analysis  for  existing  aircraft  structure  incorporating  honey¬ 
comb  sections  os  proposed  in  Phase  I  of  Reference  ($. 

The  test  program  consisted  of  two  parts.  The  first  parses  outlined  hr 
Appendix  A,  provides  test  data  for  correlation  of  a  general  buckling  equa¬ 
tion  for  honeycomb  panels  subjected  to  thermal  blast.  The  second  part, 
which  is  presented  in  Appendix  Ef,  was  conducted  to  substantiate  methods 
for  predicting  the  stress  distribution  in  a  typical  compact  te  structure  sub-  - 
jected  to  load-temperature  environment. 

The  test  methods  and  instrumentation  used  in  both  parts  of  this  test  were 
designed  to  provide  the  most  accurate  test  results  possible  based  on  exper¬ 
ience  gained  during  previous  thermal  tests  and  research.  The  most  valuable 
experience,  especially  in  test  instrumentation,  was  obtained  during  the 
thermal  tests  conducted  for  the  FJ-4B  "Hordtack  Project"  of  Reference  (S^» 
Problems  were  encountered  during  the  "Hordtack  Project"  in  devising  an 
adequate  strain  gage  and  thermocouple  system.  A  strain  gage  installation, 
which  consisted  of  nichrome  gages  and  ceramic  cement  for  readings  to 
600°F.,  proved  to  be  inadequate.  The  ceramic  cement  acted  as  an  insulator 
on  the  heated  surface  which  produced  local  stress  concentrations  due  to 
thermal  gradients.  Investigations  were  conducted  on  several  strain  gage 
systems  and  a  compromise  of  temperature  reducing  the  maximum  to  400°F. 
was  found  to  be  necessary  in  order  to  use  a  phenolic  cement  with  the 
nichrome  gage.  This  system  provided  acceptable  results  and  was  adopted 
for  future  elevated  temperature  tests.  Another  problem  was  encountered 
in  obtaining  a  positive  thermocouple  installation  on  aluminum  surfaces. 

The  previous  practice  involved  welding  the  chromel  alumel  wire  to  the 
aluminum  surface  (generally  non  weldable  7075  material).  Excessive  labor 
time  required  to  replace  the  thermocouples  necessitated  the  investigation  of 
other  methods  for  installation.  A  successful  method  which  was  adopted 
consisted  of  drilling  a  small  hole  In  the  surface,  installing  the  thermocouple 
wire,  and  peening  It  in  place.  The  400°F.  strain  gage  and  the  peened 
thermocouple  systems  were  used  in  the  test  in  this  chapter. 

The  recording  and  temperature  controlling  devices  used  during  this  are 
considered  to  be  some  of  the  best  thermal  testing  equipment  available  at 
this  time.  The  Sanborn  Recorder,  which  was  used  during  the  honeycomb 
panel  test,  is  a  direct  writing  instrument  which  can  record  a  number  of 
variables  simultaneously.  The  final  record  is  permanent  and  in  true  rectan¬ 
gular  coordinates.  The  Gilmore  Data  logger,  used  in  recording  load, 
temperature,  deflection,  and  strain  on  the  composite  structure,  is  capable 
of  printing  out  100  channels  of  test  data  in  a  total  elopsed  time  of  20 
seconds.  The  Research  Inc.,  Ignitron  Unit  used  in  controlling  temperatures 
on  both  tests,  is  a  high  quality  temperature  controller  which  electronically 
senses  and  controls  error  in  specimen  temperature. 
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The  following  test  instrumentation  and  procedure  ore  recommended  for 
combined  food  and  thermal  testing.  tn  preparing  the  specimen  for  test; 
ins  toff  Totnofl  Strain  Gages  (C12- T42  forafumtmmi)  with  GA-50  pfienoffc 
cement  for  an  adequate  strain  measuring  system  up  to  400°F.  After  In¬ 
stallation  per  manufacturers  recommended  procedures,  the  gages  ore  to  be 
cured  for  two  hours  at  350°F.  and  one-half  hour  at  450°F.  for  the  best 
results.  At  the  present  time,  strain  gages  designed  for  temperatures  higher . 
than  4Q0°f.  have  produced  unreliable  results.  The  temperature  on  the  test 
specimen  can  be  measured  or  controlled  by  chrome!  alumef  thermocouples. 

It  is  recommended  that  these  thermocouples  be  installed  at  the  strain  gages 
far  strain-temperature  correlation  and  at  locations  on  the  specimen  where 
a  certain  temperature  or  heat  rate  Is  desired.  The  thermocouples,  which 
establish  a  certain  temperature  or  heat  rate,  are  used  Irr  oanjuncflan  with 
a  temperature  controller  unit  such  as  the  aforementioned  Research  fnc. 
Ignitron  to  automatically  regulate  power  input  to  the  heating  elements. 

All  thermocouples  are  to  be  installed  on  aluminum  by  drilling  holes  for 
the  wires  and  peening  them  fast  ct  the  surface.  After  instrumenting  the 
test  specimen,  the  strain  gages  were  coated  with  a  black  Ml L-  L- 1 9537 
acrylic  lacquer  to  prevent  gege  shorts  and  the  entire  heated  area  is  sprayed 
with  an  aqua-dog  suspension  for  heat  absorbsion  and  stray  potential  shielding. 
The  quartz  lamp  radiant  heating  units  such  as  the  standard  ALT8-612  unit 
are  used  in  most  cases  to  provide  a  heating  system;  however,  for  unique 
specimen  configurations  or  in  restricted  areas,  tubular  (Cal-Rod)  beaten 
may  be  used. 
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6.0  DISCUSSION  OF  OPTIMUM  DESIGN  PROCEDURES 

The  preceding  sections  of  this  report  have  dealt  with  thearw  lysis  of  aircraft 
structures  subjected  to  temperature  gradients,  fn  the  post  it  has  been  the 
practice  to  design  an  aircraft  for  normal  flight  conditions  and  then  analyze 
the  aircraft  to  establish  the  limit  temperature  conditions  it  would  withstand. 

This  method  is  completely  unsatisfactory  in  many  cases,  specifically  for  the 
case  of  determining  the  vulnerability  of  a  weopon  system  structure.  There- 
quired  performance  of  the  weapon  system  can  only  be  designed  Into  the 
structure,  not  analyzed  into  it. 

The  work  presented  in  this  section  f£  a  first  attempt  fa  establish  an  optimiMr 
design  procedure  for  aircraft  structures  with  a  temperature  gradient.  The 
results  of  this  study  have  been  gratifying  in  that  a  relatively  simple  procedure 
has  t  een  developed.  This  section  presents  design  procedures  for  two  different 
types  of  construction  techniques.  The  first  presented,in  Section  6.1, is  for 
a  honeycomb  panel  with  fixed  load  and  fixed  strain  conditions.  These  panel 
optimization  methods  can  be  used  for  designing  the  isolated  honeycomb  panels 
for  a  given  load  or  the  non-structural  honeycomb  panel  with  a  given  strain 
induced  by  the  primary  structure.  The  more  detailed  requirements  for  these 
design  methods  will  be  presented  in  Section  6.1. 

The  two  honeycomb  panel  methods  just  outlined  are  also  used  as  a  part  of  a 
honeycomb  box  beam  design  procedure.  The  design  of  a  box  beam  to  resist 
thermal  gradients  can  be  approached  through  two  structural  concepts.  The 
first  is  to  use  conventional  construction  and  design  for  the  induced  stresses 
which  may  require  extra  material  and  therefore  extra  weight.  The  second 
method  is  to  reduce  the  thermal  stresses  as  much  as  possible  with  corrugated 
spar  webs.  The  second  method  should  be  lighter  in  weight  than  the  conven¬ 
tional  construction  but  the  induced  thermal  deflections  are  larger  than  those 
in  the  conventional  structure.  In  each  of  these  cases  the  box  beam  section 
is  optimized  for  given  applied  primary  and  secondary  moments  with  some  limit 
heat  input.  The  secondary  moment  is  the  critical  moment  in  the  reverse 
direction  from  the  primary  moment.  The  honeycomb  box  beam  optimization 
procedure  computes  the  minimum  weight,  producible  configurations  for  the 
upper  and  lower  skin  panels  as  well  us  the  spar  spacing. 

A  second  box  beam  optimization  method  was  developed  to  determine  the 
minimum  weight  configuration  for  skin  stringer  type  construction.  This  method 
includes  the  two  types  of  spar  construction  as  mentioned  in  the  discussion  of 
the  honeycomb  box  beam  method.  The  skin  stringer  box  beam  design  uses  two 
spars  and  integrally  machined  skin  stringers.  The  skin  thickness. and  stringer 
spacing  for  the  upper  and  lower  skins  ore  the  prime  objectives  in  the  optimi¬ 
zation  procedures. 

The  method  applied  in  these  optimization  procedures  have  not  used  the  allowable 
stress  methods  os  originally  planned  .  The  basic  idea  of  the  allowable  stress 

method  is  used  however,  in  that  the  design  methods  are  based  upon  the  element 
strain  rather  than  stress. 
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HTsfoncally,  optimization  methods  were  relatively  simple  when  only  a  few 
pc ra meters  were  involved,  in  many  cases  the  optimum  configuration  could  be 
found  directly  from  the  differentiation  or  minimization  of  a  weight  equation' 
which  was  expressed  in  the  terms  of  the  airowoofe  equations.  With  attempts 
to  optimize  more  complex  structural  arrangement*  ‘he  weight  equation  become 
too  complex  to  permit  direct  solution.  At  this  point  the  high  speed  computer 
became  a  tool  to  minimize  the  weight  equation  through  calculation  of  several 
configurations  and  selecting  the  minimum  configuration  from  a  simple  equation 
through  that  points.  With  the  introduction  of  temperature  gradients  info  the 
problem  the  interdependencies  of  the  structural  parameters  becomes  even  more 
complex.  It  becomes  extremely  difficult  to  ,olve  this  problem  in  the  terms  of 
applied  stress  and  allowable  stresses  os  post  optimization  problems  have  Been 
handled.  Therefore,  it  became  necessary  to  develop  a  new  procedure  or  con¬ 
cept.  If  was  known  that  by  assuming  the  total  strain  of  one  element  the  strain 
of  the  other  elements  could  be  determined  independent  of  the  element  areas 
for  an  axially  loaded  member.  Thus  several  configurations  can  be  determined 
from  different  assumed  strains  to  obtqin  the  strain  at  which  the  minimum  weight 
configuration  occurs.  More  details  of  this  technique  will  be  discussed  in  the 
following  sections  of  this  report.  With  the  addition  of  bending  info  the  pro¬ 
blem,  an  expression  for  the  bending  strain  must  also  be  included  in  the  element 
strains.  This  was  accomplished  by  assuming  a  maximum  bending  strain  on  an 
element  with  a  straight  line  distribution  between  elements  (plane  sections  re¬ 
main  plane).  This  assumed  bending  strain  is  iterated  until  a  balance  with  the 
bending  moment  is  achieved.  Again  several  configurations  are  computed  and 
a  minimum  weight  configuration  determined  using  these  points. 

The  optimum  configuration  thus  derived  may  then  be  used  as  a  basis  of  the 
final  design  development.  The  final  design  then  must  be  checked  against  the 
design  limits  using  the  allowable  stress  program  presented  in  Section  2.0  of 
this  report.  Any  changes  in  the  configuration  during  the  final  design  devel¬ 
opment  phase  may  be  referenced  to  the  non-optimum  configurations  generated 
during  the  optimization  process.  A  reference  to  these  configurations  will 
indicate  the  weight  penalties  involved  in  the  changes  and  possibly  indicate 
the  direction  in  which  the  change  should  be  mode. 
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6.1  HONEYCOMB  COMPRESStC*. '  'ANEL  OPTIMIZATION 

This  section  discusses  the  optimization  of  honeycomb  panels  (Rg.  &i  T)  at 
individual  structural  dements.  Two  different  conditions  were  amiArd 
for  the  honeycomb  panel  design.  A  method  for  the  optimum  design  of  a 
honeycomb  compression  panel  for  a  given  loco  .vas  first  derived.  A  second 
honeycomb  panel  design  method  was  derived  for  a  given  axial  strain.  The 
fixed  strain  procedure  lends  itself  to  the  design  of  secondary  structural 
elements  dependent  upon  the  deflection  and  strains  c f  adjacent  primary 
structure. 

Bath  design  methods  have  been  derived  for  thermal  gradients  dependent 
upon  the  panel  configuration.  The  dependency  of  the  temperature  dis¬ 
tribution  through  the  panel  on  facing  thickness,  core  density,  and  cor*  .. 
depth  insures  an  accurate  design  solution,  but  also  complicates  the  pro¬ 
blem  greatly.  To  reduce  this  complexity  the  facing  temperatures  have 
been  expressed  in  on  emperica!  tabular  form  as  a  function  of  the  configur¬ 
ation  parameters.  With  this  approach  a  minimum  amount  of  time  is 
expended  computing  temperatures.  However,  since  the  final  temperature 
is  dependent  on  the  configuration  and  the  configuration  dependent  on  the  ■ 
temperature  an  iterative  process  is  required  to  achieve  the  final  balanced 
design  condition. 

6.1.1  HONEYCOMB  COMPRESSION  PANEL  OPTIMIZATION  WITH 
FIXED  LOAD 


The  honeycomb  compression  pane!  with  a  fixed  load  is  designed  to  satisfy 
three  failure  modes.  The  three  failure  modes  used  in  this  derivation  are 
general  panel  buckling,  face  wrinkling,  and  interceli  buckling.  The 
genera!  buckling  allowable  load  as  used  in  the  following  derivation  is  a 
function  of  two  unequal  facing  thicknesses,  core  depth,  core  density,  edge 
fixity  conditions,  and  panel  size.  This  equation  con  be  written  as  follows: 


N  =  pM  A  H* 


(6.1.1) 


A  -  n  ^  tf  1 1  ^2  ^2 
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( figure  6-2) 


The  constants  cj,  c~,  c3,  c^  are  a  function  of  panel  edge  fixity,  a/b  ratio 
and  number  of  buckling  waves.  The  constants  can  be  computed  os  fallow*, 
for  the  four  primary  boundary  conditions: 

Simply  supported  panel 

‘2=’  *■•3- 


Loaded  edge  S.S.,  Unleaded  edge  clamped 
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The  foce  wrinkling  mode  of  fcilore  is  primarily  a  function  of  the  core  density 
and  focing  mocfufus.  The  allowable  facing- stress  is  expressed  as  follows: 

Few  -  -5  pMpLj2  GCEC  (Es  +  eE,)l,/3  (6.1.5) 

fe  exist  be  noted  at  this  point,  however,  that  the  facing,  moduli,  ^  and  Ep 
ore  a  function  of  Fcw  when  Fcw  is  above  the  proportional  limit  for  the 
esrferlal.  In  such  a  cose  the  Ramberg- Osgood  stress  strain  relation  (Section 
7.0}  must  be  used.  The  mfercell  buckling  is  <t?jo-  dependent  upon-  the 
Ramberg- Osgood  strain  relation  along  with  facing  thickness,  and  cell  size.. 
The  allowable  interceli  buckling  facing  stress  is  written  as  follows: 

3/2 

Fcj  *  .9  E^JLj  (6.1.6) 


er  =  AA!L 


E  +  E, 


(6.1.7) 


With  the  allowable  Equations  (6.1.1),  (6.1.5)  and  (6.1.6)  the  design  of  a 
honeycomb  compression  panel  can  be  accomplished.  The  requirements  of  the 
design  for  this  derivation  has  included  a  temperature  gradient  as  well  as  the 
axial  compression  load.  The  applied  stress  und^  these  conditions  can  be 
expressed  for  the  critical  compression  os  follows: 


Fc  applied  =  <*E  AT 


-  xZ 


a  E  AT  A  A  +- 


N 

£A  A 


(6.1. 


Because  of  the  two  facings  being  a{  different  temperatures  and  thus  having 
different  material  properties  and  stress  level,  an  effective  area  of  each 
element  must  be  used.  The  applied  stress  Equation  (6.1.8)  is  made  up  of 
the  full  fixity  slress  aE  AT,  the  axial  relief  stress  1/A E  .  AT 
A  A  ,  and  the  applied  load  stress  N/  £  A  A.  No  bending  relief  stress 
is  used  in  this  analysis  since  it  is  assumed  that  the  panels  are  restrained  in 
bending  by  adjacent  structure.  The  above  Equation  (6.1,8)  can  now  be 
rewritten  as  follows: 


F^opplied  =  “j  Es,  (T,-Tr)  - 
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(6.1.9) 
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For  the  design  of  the  .Minimum  weight  pane?  the  allowable  face  wrinkling 
and  general  panel  buckling  stresses  will  be  equal  to  the  applied  load 
stress.  Therefore  a  relation  for  core  density  (Ws)  can  be  derived  by 
combining  Equations  (6.  T.5)  and  (6.1.9)  as  fouawsr- 


ffl  +  Esl  ff2 


r«i  Esi  (TrV  p- 


*c  = 


|n-*2  E>2  <VV  *f2 1 

rtv) 


(6.1.10) 


ffl 


4  '*  [N 


Esl _ 

f066  Gc  Ec  (Es)  ♦  3  E„)j  1/3 

■] 


2  Es2  (T2'Tr)  *f2 


.5  [.066  Gc  Ec  (Esl  +  3  Eyl)  j  V* 


(6.1.11) 


Another  parameter  which  must  be  determined  is  the  core  depth  (C).  An 
expression  for  (C)  may  be  derived  directly  from  the  general  buckling  Equa¬ 
tion  (6.1.1)  as  follows: 

N  =  PM  A  H2 
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x  =  >H  *  *f2 
2 


Y  =  ffi  ff2  Etl  ^2 

ffl  Etl  ^  ff2  ^2 


Substituting  Equations  (6.1.2).  ((5.1.3)  and  (6.1.4)  >  N  become*: 

A  +P5i-  +  V  ^  F  \ 

U  V c4.  yi  \  KjY  (C  +  XT 


/K+lt-  +  M  F\ 


(6.1.12! 


KjYC 


vx  “  c# 
cx 


K|Y  C 


G»  =  G  *y  for  square  cell 

CX  '  (Figure  6.2) 

Got  *  Gc//2  for  Hexcell  . 

KKj  (C  +  X)2  y/KjYC  +  KjYCN  F^YfC+X) 


K,2  Y2  C2 

1  +  L  + - 

Ci  G*  GjT 
4  cx  cy 

K1  K  (C  +  X)2  +  C  (C  +  X)2/Kl2 


N_  +  J^L  +  C2  /nk^fy 
Y  Y 

+  (C3  +  2C2X  +  C  X2)  [7^ 


C  KiK  +  2CXKjK  +  X2K,K 
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Alf  the  Equations  (6.1.  II),  (6.  1.6)  (6 .1.17)  for  the  design  of  a  honeycomb 
panel  in  compression  have  now  been  presented.  From  these  equation*  it  It 
evident  that  the  design  procedure- for  a  given  panel  with  a  thermal  input  Is  vny 
complex,  if  attached  directly.  The  complexity  is  primarily  established  by  the 
Infer  ^dependency  of  the  temperature  distribution  with  the  panel  configurations 
and  the  inelastic  stress-strain  relation  inherent  in  an  efficient  design. 

While  designing  a  honeycomb  panel  to  perform  a  certain  {ob  it  will  be  found 
that  there  ore  on  infinite  number  of  panel  configurations  which  will  perform 
this  Job  with  a  zero  margin.  The  structural  efficiencies  for  these  infinite 
number  of  configurations  vary,  however,  over  a  very  wide  range.  The  con¬ 
figuration  desired  is  the  one  which  attains  the  best  efficiency  (or  minimum 
weight). 

The  following  design  procedure  has  been  derived  to  compute  this  minimum 
weight  (optimum)  configuration  under  the  conditions  previously  stipulated. 

The  basic  approach  to  this  problem  has  been  changed  from  one  of  stress  as 
used  in  the  past  to  one  of  strain.  By  assuming  the  total  strain  on  one  facing 
of  the  panel  the  strain  can  be  computed  for  the  other  facing.  The  strain  on 
the  second  facing  is  expressed  as  follows: 

«2=  *1  -  •  (TrV 


This  equation  Is  valid  since  it  has  been  assumed  that  no  bending  occurs  in 
the  panel.  The  bending  restraint  is  assumed  to  be  provided  by  adjacent 
structure. 


At  this  point  an  arbitrary  facing  thickness  tfj  is  also  assumed.  With  this 
assumed  value  of  tf  j  an  initial  estimate  of  the  facing  temperature  is  com¬ 
puted.  This  Initial  estimate  is  made  with  a  very  simple  equation; 


Tj  =  +  ^2  r  =  Boundary  layer  Temperature  (6.1,18) 


Knowing  and  82  and  the  temperature  of  the  facings,  the  facing  stress 
OO  can  be  computed  using  the  Romberg- Osgood  equation.  The  solution  of 
the  Ramberg-Osgood  equation  is  an  iterative  process  of  trial  and  error.  The 
equation  to  be  solved  is  as  follows: 


(6.1.19) 


The  yield  stress  (Fy)  is  obtained  from  material  properties  at  temperature. 

Knowing  the  tjj  end  the  facing  stresses  Fcj  one!  Fc2  the  second  facing  thick¬ 
ness  can  be  computed  os  a  function  of  the  applied  load  (N).  The  expression 
for  tf2  can  be  written  os  follows: 
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*rz~  K ~ *nFc? 

fc2 


(6-1.20) 


Having  obtained  (Equation  6.1.20)  two  other  parameters  remain  to  be 
determined  to  complete  the  panel  configuration.  These  parameters  tc/s, 
and  C,  can  be  determined  from  Equations  (6.1.11)  and  (6.1.17),  respectively. 
With  the  final  determination  of  panel  configuration  for  the  given  lood,  temp¬ 
erature,  tfj,  and  8]  a  new  estimate  of  the  temperature  distribution  can  be 
made.  This  new  estimate  can  be  made  with  much  greater  occurocy  because 
a  panel  configuration  has  been  determined.  The  method  of  determining  this 
temperature  distribution  is  taken  from  Reference  10.  The  curves  used  for 
aluminum  honeycomb  are  presented  in  Figures  6.3  and  6.4. 

The  process  of  estimating  the  temperature  and  recomputing  the  panel  con¬ 
figuration,  as  described,  must  be  repeated  until  temperature  convergence 
Is  achieved.  At  this  time  the  temperature  distribution  and  the  panel  con¬ 
figuration  have  become  compatible.  The  area  cr  relative  weight  of  this 
configuration  can  be  determined  as  follows: 

A=  ifl  +  *f2  +  *  tc  C  (6.1.21) 

S 


where 


K  =  2  for  square  cell,  3  for  Hexcell 

With  a  compatible  panel  configuration  determined  for  the  assumed  strain 
8j  and  facing  thickness  to  another  strain  8  will  be  assumed.  A  com¬ 
patible  pcne!  configuration  is  then  determined  for  this  assumed  strain  also. 
This  process  is  repeated  for  a  third  time  so  that  there  are  three  different 
configurations  for  the  given  load  and  temperature  condition  and  assumed 
facing  thickness  tfj.  At  this  time  the  strain  at  which  the  minimum  weight 
configuration  occurs  can  be  determined.  The  process  of  determining  this 
minimum  is  one  of  fitting  an  equation  (described  in  Section  6.4)  for  area  as 
a  function  of  5.  through  the  three  computed  points  and  differentiating  the 
equation  and  setting  it  equal  to  zero.  The  5j  which  satisfies  the  dA/d  8j  = 
O  equation  is  the  strain  at  which  the  minimum  weight  (or  area)  configuration 
occurs.  Using  this  8j  the  optimum  configuration  for  the  assumed  t^j  can  be 
computed  by  the  preceding  method. 

At  this  point  in  the  analysis  the  assumed  tj  j  is  changed  and  the  complete 
procedure  of  computing  three  panel  configurations  for  three  different  5  j's 
is  .epeated.  Again  following  the  third  5^  computation  the  optimum  con¬ 
figuration  is  determined  for  the  second  assumed  tfj.  This  process  is  repeated 
for  a  third  tn  to  establish  three  optin.jm  8  j  configurations  as  o  function  of 
tf]  (Figure  6.5) 
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With  these  three  points  the  absolute  optimum  con  be  determined.  The 
determination  of  the  absolute  optimum  point  cut  be  accomplished  again  by 
fitting  an  equation  for  the  area  (described  in  Section  6.4}  through  the  three 
computed  points  as  a  function  of  tfj  and  setting  the  derivative  (dA/dtf|)  equal 
to  zero.  After  solving  this  equation  for  the  optimum  tfj  the  process  must  be 
repeated  to  determine  the  configuration  compatible  with  the  optimum  tf]. 

This  final  configuration  will  satisfy  the  loading  and  thermal  conditions, 
applied  with  the  minimum  weight  configuration,  within  the  toferancekeC 
engineering  accuracy. 

The  preceding  method  was  programmed  as  described  on  the  IBM  709  Digital 
Computer  as  Program  E2001 .  Details  of  this  program  such  at  the  FORTRAN 
fisting,  block  diagram,  symbol  listing  ,  end  input  data  ore  presented  in 
Appendix  D  . 

Results  from  this  program  are  plotted  in  Figure  6.6.  Each  configuration 
computed  in  the  described  procedure  is  plotted  to  present  a  better  picture 
the  process  to  find  the  optimum  configuration  that  satisfies  the  given  loading 
and  heating  condition.  This  plot  is  also  a  valuable  tool  in  selecting  non- 
optimum  configurations  in  a  case  where  the  optimum  configuration  cannot 
be  used. 

6.1.2  HONEYCOMB  PANEL  OPTIMIZATION  WITH  FIXED  STRAIN 

The  fixed  strain  design  procedure  has  been  derived  to  satisfy  the  requirements 
of  secondary  structural  elements.  These  secondary  structural  elements  in 
general  are  not  highly  stressed  lead  carrying  elements  but  elements  which 
ore  required  to  carry  only  the  strain  induced  in  them  by  the  primary 
structure  and  thermal  loadings.  However,  these  secondary  structural 
elements  must  still  be  capable  of  withstanding  the  induced  stresses  of  strain. 
The  failure  modes  which  establish  the  design  criteria  are  general  panel 
buckling,  face  wrinkling,  and  intercell  buckling.  The  Equations  written 
for  these  failure  modes  in  Section  6.1.1  remain  basically  the  same,  (6.1.1) 
(6.1.5),  (6.1.6).  The  approach  to  the  derivation  of  this  method  remains 
one  of  strain  analysis.  If  5  jg  is  considered  as  the  +  extensiona!  strain 
induced  by  the  primary  structure  the  total  strain  for  each  of  the  two  facings 
can  be  written  as  follows: 

*i  - »,  <vy  +  SN 

8 2  ~  S  j  *  (°iT|”°2^y  (6.1.23' 

The  initial  temperature  estimate  is  made  by  the  following  equation 

T|  =  T2  +  T2  =  Tb  (6.1.24) 

»fl 
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Knowing  5  ^  and  5~  the  second  modulus  and  element  stress  for  each  facing 
eon  be  found  from  the  rtx 


:  material  stress  strain  curve  (Reference  Section- 2.0). 
At  this  point  the  core  configuration  can  be  computed] 
facing  thickness  tf  j  with  the  following  equation: 


*c  r 

—  =*  _ 

s 

;  >  •< 
1 


S'  1^1 


s*5- 

s  =  Af_£Ri_\  ^ 

\Ki  /  »n 


V* 


3 


3/2 


(£.1.26) 


The  second  facing  thickness  tf2  can  be  computed  from  the  equation  for  cell 
size  (S)  as  follows: 


tf2  =  S 


(^— ) 

\-9Er2  / 


2/3 


(6.1.27) 


The  remaining  panel  parameter  to  be  computed  is  core  depth  which  can  be 
derived  from  the  general  buckling  equation: 


N  = 


K,  (C  +  X’,[K  ^ 


_S|xC4_ 
_1_  +  L  +  K}2C2 
Y  G^G»yc4 


-(6.1.28) 


The  load  N  to  be  sustained  by  the  panel  can  be  expressed  os  follows: 
NA=  Esl  s\  *fl  +  Es2  *2  *f2  (6.1.29) 


Setting  the  allowable  load  (Equation  6.1 .28)  equal  to  the  applied  load 
(Equation  6.1.29)  the  following  cubic  equation  for  the  core  depth  C  is 
obtained; 


=  _L  x2k,k 

Y  NA 


(Cl  + 


1  - 

~r~ 


C2) 


(6.1.30) 


*1 _ + 

Gcx*4 


(C3  + 


1  -  F 


C2  a 


(6.1.31) 


This  cubic  equation  for  C  is  solved  by  a  trial  and  error  iteration  method  on 
the  computer.  At  this  point  in  the  computation  the  temperature  distribution 
in  the  pane!  must  be  computed  and  compared  for  convergence  with  the  last 
computed  value.  The  temperature  information  is  taken  from  the  curvesin 
Figures  6.3  and  6.4.  If  the  temperature  has  not  converged  the  panel  con¬ 
figuration  must  be  recomputed  as  described.  Following  convergence  of  the 
temperature,  the  area,  or  relative  weight,  of  the  panel  can  be  computed 
from  the  following  equation: 

A  =  tfi  +  tf2  +  K  i^L  (6.1.32) 

where  K  =  2  for  square  cell,  3  for  Hexcell  (See  Figure  6.2) 

The  process  described  is  repeated  for  several  values  of  tfj  for  the  given 
applied  load  strain  8n.  This  procedure  produces  a  series  of  points  re¬ 
lating  the  panel  area  (or  weight)  as  a  function  of  tfj.  It  was  found  that 
this  procedure  did  not  produce  a  minimum,  but  went  to  zero  Area  for 
tf|  -  O.  At  this  point  in  the  study  minimum  values  for  the  facing  thicknesses 
tfi  and  tf2  were  used  along  with  core  density  and  cell  size  limits.  With 
these  limits  imposed  a  minimum  practical  configuration  is  reached. 
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The  method  discussed  wor  programmed  an  the  IBM  709  Digital  Computer  * 
at  Program  E2004.  Details  of  this  program  such  as  the  FQ2I2AN  1'itiep. 
block  diagram,  symbol  listing  and  input  data  are  presented  in  Appendix  D. 

Sample  results  from  this  program  are  plotted  in  Figure  6.7.  Each  configura¬ 
tion  computed  in  the  described  procedure  is  plotted  to  present  a  better 
picture  of  described  process.  Each  configuration  derived  will  satisfy  the  - 
induced  strain  and  temperature  requirements  of  the  design  condition. 
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6.2  HONEYCOMB  BOX  BEAM  OPTIMIZATION 

To  this  point  the  optimization  of  honeycomb  structural  elements  hove  been  • 
discussed,  hr  general  rfresy elements  appear  m  a  primary  food  carrying 
structure  which  can  be  classified  as  a  box  beam.  The  honeycomb  elements 
discussed  to  this  point  ore  used  in  this  box  beam  as  the  skins  or  facings  of 
the  box.  The  spor  webs  and  spar  spacing  must  be  designed  to  satisfy  the  design 
condition  as  well  as  the  facings.  While  the  honeycomb  panel  optimization 
design  condition  allowed  for  only  axial  hading,  the  box  beam  must  he. 
designed  for  bending  and  shear  loads.  The  bending  load  can  stiff  be  assumed 
to  induce  axial  loads  in  both  the  upper  and  lower  honeycomb  facing  panels. 
Besides  the  applied  bending. moment  there  is  now  a  thermally  induced  bending 
moment  to  be  accounted  for.  This  thermal  moment  is  dependentupon  the 
upper  pane!  configuration,  which  determines  the  primary  tempe  rafure,  and  the 
remaining  portion  of  the  box  beam. 

The  optimum  design  process  developed  during  this  study  uses  both  the  fixed 
load  and  fixed  strain  honeycomb  panel  procedures  to  derive  the  optimum  box 
beam  configuration.  The  details  of  application  will  be  discussed  later  in 
this  section. 

In  the  design  of  a  box  beam  two  types  of  spar  web  construction  must  be  con¬ 
sidered  when  the  thermal  stress  problem  is  present.  These  two  web  configur¬ 
ations  are  designated  as  straight  web,  Figure  6.8  and  corrugated  web. 

Figure  6.9.  The  two  types  of  construction  differ  drastically  in  the  presence 
of  thermal  stresses.  In  the  case  of  the  straight  web  box  beam,  thermal 
expansions  of  the  upper  facing  induce  loads  in  the  spar  v/eb  and  lower  facing 
as  exicnsiona!  forces  and  induced  bending.  In  the  case  of  the  corrugated  spar 
web  the  upper  facing  panel(including  the  spar  attachment  channels  ) Is  Free 
to  expand  without  inducing  loads  into  the  rest  of  the  structure.  Thus  the 
thermal  effects  are  retained  entirely  within  the  hot  facing  panel.  There  are 
certain  advantages  to  this  type  of  structural  arrangement  but  the  major  dis¬ 
advantage  is  the  large  deflection  induced.  If  this  increased  deflection  and 
warping  under  thermal  loads  can  be  tolerated  a  light  weight  structure  is 
most  likely  possible. 
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CORRUGATED  WEB  BOX  BEAM 
FIGURE  6.9  . 

6.2.1  SPAR  WEB  DESIGN 


In  the  present  study  two  types  of  spar  webs  are  used  for  the  box  beam  design 
The  first  type  to  be  discussed  is  the  straight  spar  web  design  to  carry  the 
average  axial  strain  and  shear  strains  without  buckling.  The  allowable 
shear  stress  can  be  expressed  as  follows: 


1  scr 


*  Kws 


/»w\  2  p 

Et* 


(6.2.1) 


The  allowable  axial  stress  is  expressed  as: 
tw  .  2 


^cr 


•w 


(6.2.2) 


The  applied  sheai  stress  and  axial  stresses  are: 

=  Q _ 

fs  /.£k  +  l] 

Iw 

/CL 


1]  ^bw  w  ew^w 

+-  =  No.  of  Spar  webs 


(6.2.3) 


These  values  are  combined  in  the  shear  and  axial  load  combined  stress 
equations  to  derive  an  expression  for  t  as  follows: 


5T? 
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Solving  for  tw  : 


This  equation  (6.2.6)  can  now  be  put  in  the  form  of  a  simple  cubic  if  t  £ 
is  set  equal  to  X^,  which  can  be  solved  directly  for  Xfy,  or  t  £. 

The  second  type  of  spar  web  design  is  the  corrugated  web.  The  type  of 
corrugation  selected  for  this  study  is  shown  in  Figure  6.10. 


CORRUGATED  SPAR  WEB 
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Sinc-i  this  spar  webis  incapable  ctf  carrying  any  part  of  the  bending  moment 
on  the  box  beam,  only  the  shear  food  most  be  considered  in  Its  design.  This 
type  of  shear  web  can  be  designed  to  an  ultimate  shear  vafue  (FjyJ  selected 
by  the  designer.  The  corrugation  length  (bc)  can  then  be  determined  to 
prevent  buckling  of  the  web.  The  web  design  parometeis  con  be  computed 
as  follows: 


a 


CL  +  1  \ 

w  ) 

If  '<  min 
bc  min  =  minimum  rivet  spacing 


Bearing  Load  -  tw  BRj 
in. 


Rivet  load  = 

Rivet  Dia.  Req.  =  -^c-— — 


/KwsJw\  1/2  (6.2.7) 

bc  “  X  J 

set  bp  —  b^  min 

(6.2.8) 

(6.2.9) 

(6.2.10) 
(6.2.11) 


If  Rivet  Dia.  >  Max.  Allow  Rivet  Dia. 

Set  RD  =  MRD  b  =  RP  **"*"***  (6.2. 12) 

Q 

IF  bc  <  bc  min  tw  =  tw  bc  min  (6.2.13) 

bc 

set  bc  =  bc  min 
=  2  x  BR  x  tw 

The  relative  area  of  the  web  can  then  be  computed  as: 


(6.2.14) 
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6J2JL  STRAIGHT  SPAR  WEB  BEAM  DE51GN 

The  method  selected  for  (he  design  of  the  upper  and.  lower  cap  panels* -  - 
dependent  upon  (he- type  of  spar  weh  consfruefTon.  The  primary  difference 
in  (he  (wo  spar  web  types  considered  in  this  study  is  the  bending  effectiveness 
of  the  spar  webs.  The  straight  web  Is  co  nsidered  to  carry  a  share  of  the 
bending  moment  while  the  corrugated  web  carries  no  bending. 

With  the  spar  web  carrying  part  of  the  bending  moment  the  portion  of  the 
load  carried  by  the  cops  is  not  dTrectfy  obtainable.  Since  the  configuration 
of  the  box  beam  is  the  object  of  the  design  problem  no  moment  of  inertia  or 
s«?e  parameters,  are  available  to  estimate  the  stress-  distribution.  Tins  ts 
specifically  true  in  the  case  of  inelastic  stress  where  the  distribution  is  non¬ 
linear.  Therefore  strain  analysis,which  is  more  independent  of  the  con¬ 
figuration/must  be  used.  In  the  case  of  the  honeycomb  panels  a  strain  on  one 
face  could  be  assumed  and  the  strain  on  the  other  face  computed  as  a  function 
of  the  temperature  distribution.  This  is  no  longer  true  since  there  is  bending 
present.  However,  if  the  assumed  upper  panel  strain  is  a  total  strain  including 
even  the  bending,  a  bending  strain  distribution  for  the  remaining  elements  can 
be  assumed.  This  assumed  bending  strain  will  take  the  form  of  a  straight  line 
distribution  of  zero  on  the  upper  cap  and  a  tension  strain  elb^  on  the  lower 
cap.  At  this  time  the  total  strain  on  each  element  can  be  computed  as 
follows: 

Average  spar  web  strain 

«w  =  eul  “  <«  0  Tu1  -  «  w  Tw)  +  1/2  (6.2.16) 

Lower  cap  strain 

eL  =  eu1  "{®o  Tul  "alTl?  +  eU>  (6.2.1?) 

Parameters  which  must  be  given  for  the  solution  of  this  problem  are  as  follows: 


CL 

Chord  length 

h  or  bw 

Beam  depth  to  center  of  faces  ' 

Primary  bending  moment 

Mr 

Secondary  reverse  bending  moment 

Q 

Primary  shear  load 

L 

Rib  spacing 
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The  initially  assumed  parameters  ate  as  follows: 

*yj  Total  strain  on  of  oppsr  aep 

ejj,  Bending  strain  on  lower  cap 

*|j  Upper  facing  thickness  for  upper  panel 

W  Spar  spocnvj 

The  temperature  of  the  upper  panel  face  (Tul)  used  in  the  element  strain 
equations  Just  described  was  token  from  the  fixed  strain  honeycomb  panel 
optimization  procedure  described  in  Section  6.1.2.  The  area  of  the  lower 
facing  is  estimated  by  the  following  equation: 


+  tFS  = 


(6.2.18) 


To  complete  the  Initial  sizing  of  the  configuration  a  spar  web  thickness  (♦w) 
must  be  obtained.  The  sixth  root  polynomial  derived  in  Section  6.2.1 
(Equation  6.2.6)  can  now  be  solved  as  a  function  of  e_. 


4  “  *w  <w  hi  -j 
Kl 


i  ~[  Q  byy 

1  V(£-  +  ,J 


(6.2.19) 


Equation  (6.2.19)  designs  the  straight  spar  web  to  remain  unbuckled  under 
the  shear  and  bending  and  axial  stresses  induced  by  the  thermal  and  primary 
static  loadings.  Now  that  a  configuration  has  been  derived  the  equili¬ 
brium  of  the  section  must  be  checked.  The  first  requirement  for  this  equili¬ 
brium  check  is  to  determine  the  location  of  the  neutral  axis,  since  •  . 
the  moment  must  be  computed  about  the  neutral  axis  of  the  beam  cross- 
section.  The  location  of  Hie  neutral  axis  below  the  center  line  of  the  upper 
panel  Is  confuted  horn  the  following  equation: 


2 


(6.2.20) 


-  ,  ^  (1  +1-)+  bwEiL  (tF4  »  *F5)  <6.2 

(  W  +  Cfj^^w  +  EsL  (*F4  +  *F5)  *  EsFl  If  1  +  EsF2  fF2 


(6.2.21) 
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Using  a  moment  bctocrce  cqooftort  about  the  nevtmf  exfir,  otr  ocftraT 
con  be  computed* 

eLb  =  '  Mu  +  +eu2Es2tF2CLy  + 

K\~  1E>  +  wTw*  ^-9:  •' 

K>-  >'i+  E^ijgcu;  » 

[-•5  [w,  ' 


(6.Z2Q 


F3CL(y  -  M 


With  this  computed  eit  the  upper  panel  load  required  for  axial  balance 
can  be  computed  as  follows: 


Pul  *F2  eue  Es2  CL  +  e1sl  Esl  fFl  CL_ 

"  +  Tw  bw  ew  Esw  -  e^  Es|_  *ps  CL  (6.2.23) 

The  load  P  j  is  then  used  in  the  fixed  load  honeycomb  pane!  method  to 
determine  the  optimum  upper  pane!  for  the  given  load  ond  optimum  strain 
level  eu]  for  that  load. 


Since  the  Pyj  and  have  been  computed  from  the  equilibrium  equations 
and  are  compatible  with  e^,  convergence  of  ey^  is  all  that  is 
necessary.  Therefore,  by  comparing  the  present  value  of  ey,  with  the  last 
value  convergence  can  be  determined  when  they  become  equal  within  a 
preset  limit. 


The  only  portion  of  the  total  configuration  which  remains  undetermined  at 
this  point  is  the  core  configuration  for  the  lower  box  cap.  The  area  of  the 
lower  panel  facings  has  previously  been  determined  from  requirements 
imposed  by  the  primary  bending  putting  the  lower  cap  in  tei  sion.  The  core 
configuration  will  be  determined  from  the  compression  loading  imposed  by 
the  secondary  bending  moment  and  thermal  input  on  the  lower  side.  The 
compression  load  (J'Prj)  imposed  on  the  lower  face  is  estimated  as  follows: 
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*ttl  =  XPul  =  ffFI  Fa+t^Fc^  01  (6.2.24X6.2.253 


_ '  Mr  xpul 


bw  Pu 


(6.2.26) 


To  determine  the  flrrcf  fower  configuration  three  thickness  ratios  of  ffjAf4  are 
used  knowing  ;he  value  of  t^  +  ^4  (Equation  6.2. 18).  By  assuming  an  initial 
core  density  and  deptbr  the  temperatures  T4  and  7^  can  be  computed  from  tW 
curves  of  Reference  (IQ.  Assuming  a  strain  e$,  can  be  computed  as: 


e4  “  e5  "  (  a5  T5  "  °4  Td) 

A  check  of  the  equilibrium  equation  must  be  performed: 


Pw  =  e,E. 


4£s4tF4  +  e5  Ej5  fF5 


(6.2.27) 


(6.2.28) 


the  computed  Pjq  must  be  equal  to  the  XPjg  applied  to  the  panel.  If.this  is 
not  satisfied  ej  is  changed  until  the  condition  is  met.  When  the  condition 
is  met  the  core  density  can  be  computed  as: 


•066  Grf  Erf  (Erf  +  3ETs)  J1/3 


J  (6.2.29) 


The  core  depth  (C)  can  also  be  computed  from  Equation  (6.1.17)  of  Section 
6.1.1.  With  the  computed  configuration  just  established  a  new  temperature 
distribution  in  the  panel  must  be  computed  and  checked  for  convergence, 
against  the  last  temperature  values  used.  If  convergence  has  not  been  ob¬ 
tained  the  configuration  must  be  recomputed  as  described.  When  conver¬ 
gence  is  obtained,  the  area  of  the  lower  panel  Is  computed  as: 


A5='F4+*F5+K^  C 


(6.2.30) 


The  procedure  must  then  be  repeated  for  another  ratio  of  tfs/tf4  until  three 
configurations  have  been  obtained.  At  this  time  a  minimum  weight  config¬ 
uration  can  be  estimated  with  the  curve  matching  process  (Section  6.4)  from 
the  values  A^  as  a  function  of  *f5* 

When  the  lower' panel  configuration  has  been  determined  as  discussed  above, 
the  relative  area  of  the  total  box  beam  can  be  computed  by  combining 
Equations  (6.1.21)  (6.2.30)  and  the  spar  web  area 
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CL 

W 


+T 


CL 
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(6.2.31) 


Following  the  computation  of  the  total  area  Ay  a  new  box  beam  configuration 
must  be  computed  for  another  value  of  tf  j  until  three  configurations  are 
obtained.  These  three  configurations  are  then  used  to  determine  the  value  of 
tfj  which  will  produce  the  minimum  weight  configuration  for  the  given  spar 
spacing  m - 

Following  the  determination  of  the  minimum  configuration  for  the  one  spar 
spacing  (W)  the  process  Is  repeated  for  two  more  values,  of  W  to  achieve,  a 
minimum  weight  configuration  with  respect  to  spar  spacing  (W).  This  total 
process  then  provides  the  designer  with  the  minimum  weight  box  beam 
structure  to  withstand  all  the  induced  thermal  conditions  and  loads  for  the 
many  parameters  involved. 

The  details  of  the  programmed  procedure  are  presented  in  Appendix  D  in 
the  form  of  block  diagrams,  FORTRAN  listings,  and  sample  data  sheets.  A 
sample  plot  of  data  fora  typical  design  problem  is  presented  in  Figure  6.11. 
This  plot  presents  each  point  computed  in  the  optimization  process  described 
and  presents  a  better  picture  of  the  procedures  discussed.  The  plot  also  shows 
the  weight  savings  which  can  be  achieved  by  such  on  optimization  process. 

6.2.3  CORRUGATED  SPAR  WEB  BEAM  DESIGN 


The  corrugated  spar  web  beam  design  must  be  considered  when  thermal 
gradients  are  present  in  the  structural  design  requirements.  Since  the  corru¬ 
gated  spar  web  will  not  restrain  the  hot  facing  panel  from  expanding  the  thermal 
stresses  induced  in  the  structure  are  much  lower.  These  lower  stresses  will 
permit  a  lighter  weight  design  when  considering  a  strength  design.  However, 
the  corrugated  spar  will  allow  greater  deflections  for  the  same  relative  weight, 
and  would  be  considerably  heavier  than  the  straight  web  construction  in  the 
case  of  a  stiffness  design. 

The  load  carried  by  the  upper  and  lower  caps  can  be  determined  directly 
for  the  corrugated  spar  case.  These  loads  are  computed  os  follows: 

P„!  -  P«2 

K 

This  load  can  then  be  used  directly  in  the  fixed  load  honeycomb  panel 
optimization  procedure  described  in  Section  6.1.1.  Following  the  deter¬ 
mination  of  the  upper  cap  panel  configuration  the  area  of  the  lower  cop  can 
be  determined.  The  lower  cap  area  is  assumed  to  be  critical  in  tension  for 
the  primary  bending  moment  My. 


(6.2.32) 
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FIGURE  6^11 
OPTIMIZATION  CURVE 

RIB  SPAONG-L  =  10  IN 
BOX  DEPTH- 1^  =  5  IN 
BOX  CHORD -CL  =  40  IN 
PRIMARY  BENDING  -  Mj  =  1  x  106  IN-LB 
-  SECONDAEY  BENDING- M^xlO5  IN-LB 
MAX  SHEAR  LOAD-Qj  =  100Q  LI 
RADIANT  EXPOSURE-Qj.  =  10  CAl/crn2' 
TIME  TO  PEAK,  -  n  =  1 .014  SEC 
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OPTIMUM  CONFIGURATION 
UPPER  PANEL 


tj.,1  =  .0436 

tp2  =  .0801 

C]  =  .33 


.00583 


TEMP  1  =  241°F 
TEMP  2  =  158°F 

Sj  =  .5676 


LOWER  PANEL 


tF5  =  .0499  TEMP  5  =  226°F 
tp4  =  .0334  TEMP  4  =  203°F 
C2=  .18 

-k-  =  .00326 


S2  =  .7633 
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(6.2.33) 
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Hie  lower  ponel  must  also  withstand  the  compression  load  Induced  by  the 
reverse  moment  Mg. 

The  Pg2  con  then  be  used  In  the  lower  panel  design  procedure  described 
In  Section  6.2.2,  This  design  method  ensures  sufficient  area  ut  the  kwe 
panel  to  carry  the  tension  loads  induced  bv  and  panel  stability-  ter  carry 
the  compression  loads  Induced  by  Mg. 

The  design  of  tie  corrugated  web  Is  the  only  problem  remaining  In  the  bat 
beam  design.  The  design  procedure  for  the  corrugated  web  is  described  in 
Section  6.2.1. 

The  design  procedure  to  this  point  has  been  for  a  given  tf  j  and  W.  As  In 
the  straight  web  box  beam,  two  more  values  of  tf  j  ore  used  for  the  fixed  W. 
The  area  of  each  configuration  is  computed  as; 

aT  =  fF1  +  *F2  +  fF4  +  fF5  +  K  —  Cj  +  K  ]c5  C5 


The  parabola  method  of  Section  6.4  Is  again  used  to  estimate  the  minimum 
area  configuration  as  a  function  of  tf|.  This  process  Is  repeated  for  two 
more  values  of  W  to  provide  data  to  estimate  the  minimum  configuration 
with  respect  to  tfj  and  W. 

The  IBM  program  for  the  conugated  spar  web  box  beam  design  is  included 
in  the  same  program  with  the  straight  web  box  beam  presented  in  Appendix 
D  .  Using  the  same  computer  program,  the  option  of  straight  web  or  cor¬ 
rugated  web  or  both  box  beam  designs  is  up  to  the  engineer.  Data  from  the 
corrugated  box  beam,  program  is  presented  in  Figure  6.12r  for  comparison 
with  an  equivalent  straight  web  box  beam. 
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FIGURE  6. 12 
OPTIMIZATION  CURVE 
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OPTIMUM  CONFIGURATION 
UPPER  PANEL 
»F1  =  .0411 
ff 2=  .0867 
Ci  =  .90 
»c/S  =  .  00543 


TEMP  I  =  256°F 
TEMP  2  =  119°F 


tp5  =  .0445 
Ip4  =  .0389 
C2  =  .21 
!</S  =  .002 


Sj  =  .4925 
LOWER  PANEL 

TEMP  5  =  239°F 
TEMP  4  =  215°F 


SPAR  WEBt  =  .08 


S2  =  2.5 
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63  SOX  BEAM  WITH  SKIN  AND  STRINGERS  ' 

The  purpose  of  this  section  is  to  determine  an  optimum  box  beam  structure, 
under  a  non-uniform  temperature  distribution  end  given  overall  dimensions 
and  for  a  given  applied  bending  moment,  shear  load,  and  reverse  applied 
bending  moment. 

Initially,  for  abetter  understanding  of  the  problem,  an  I-Beam  was  investi¬ 
gated.  to  obtain  same  preliminary  optimization  equations.  It  was  found  that 
fbe  best  approach  was  through  the  use  of  strain  equations.  To  start  with,  an 
upper  cap  thickness  was  assumed  along  with  an  applied  strain  and  appropriate 
strain  and  moment  equilibrium  equations  ’•'•■*tfeir.  The  strain  was  then  varied 
until  the  equilibrium  equations  were  sclir*’:  .v.  The  only  restrictions  placed 
on  the  cross-section  were  that  the  upp-  ■  c.ap  and  web  be  at  their  critical 
buckling  strains.  The  lower  cop  area  was  determined  !>/  assuming  it  is 
critical  in  tension  and  calculating  the  area  from  the  applied  moment. 

Another  thickness  was  assumed  and  the  applied  strain  varied  until  the  equil¬ 
ibrium  equations  were  satisfied.  After  several  values  of  the  upper  cap 
thickness  were  assumed,  a  plot  was  made  of  total  area  versus  upper  cap 
thickness.  The  upper  cap  thickness  at  the  minimum  area  was  then  taken  as 
the  optimum  thickness. 

It  is  this  some  general  approach  of  using  strain  equations  along  with  appro¬ 
priate  equilibrium  equations  that  is- used  in  determining  an  optimum  box  beam 
with  skin  and  stringers.  Two  types  of  box  beams  ore  considered  here,  one 
with  straight  webs  and  one  with  corrugated  webs.  The  box  beam  with 
straight  webs  will  be  considered  first. 

6.3.1  BOX  BEAM  WITH  STRAIGHT  WEBS 

If  its  assumed  that  the  optimum  structure  will  be  similar  to  the  one  shown 
beibw  in  Figure  6. 13. 


BOX  BEAM  WITH  STRAIGHT  WEBS 
FIGURE  6.13 


119 


NORTH  AMERICAN  AVIATION  ,  INC. 

ceittMtus  division 


The  structure  Ts  onafyzed  In  three  ports ;  determination  of  lower  cap  area 
and  skin-stringer  configuration  (Section  6.3. 1.1),  determination  of  upper  • 
cop  area  and  sk?n-str?nger  configuration  (Section  6J3.1.2}.  and  analysis 
of  the  box  beam  (Section  6.3.  T.3). 

6.3. 1.1  DETERMINATION  OF  LOWER  CAP  AREA  AND  SKIN-STRINGER 
CONFIGURATION 

The  lower  cap  is  acted  an  by  art  applied  moment.  My.  and  a  reverse  ~nnmrwt, 
Mr.  The  cap  has  a  total  width  CL.  It  is  assumed  that  the  lower  cqp  ts  at  a 
uniform  fenperafure  and  colder  than  the  upper  cap. 

As  a  first  approximation  it  is  assumed  that  the  lower  cap  is  critical  tn  tension. 
The  area.  A^r  of  the  lower  cap  is  then  calculated  by: 

AL  =  (6.3.1) 

kwfyuL 

where  Fj^j  j_  is  the  ultimate  tension  stress  at  room  temperature  and  bw  is  the 
arm  through  which  Aj_  acts  (in  this  case  the  depth  of  the  box  beam). 

A  minimum  skin  thickness  of  .020  inch  is  assumed  as  a  practical  limit.  There¬ 
fore,  Almost  be  >  .020  CL.  If  A|_,  calculated  by  Equation  (6.3.1)  is  less 
than  this  value.  is  taken  as: 

Al  =  .020  CL  (6.3.2J 

Since  there  Is  a  reverse  applied  moment,  Mr,  present  that  will  induce  com¬ 
pression  into  the  lower  cap,  the  lower  cop  must  be  checked  to  see  if  it  is 
critical  in  compression  and  add  the  necessary  stringers  for  the  given  area. 

A[,  to  sustain  the  compression  load.  An  effective  area,  Acff,  (that  port 
of  the  lower  cap  area  that  remains  stable  under  a  compression  load)  is 
determined  by: 

Aeff  =  Mr  (6.3.3) 

Kv  fallow 


where  Fallow  is  some  allowable  compression  stress  (Fa||ow  =  ^cy  ** 
suggested. 


The  value  of  Aeff  is  then  checked  against  A[_  calculated  by  Equ  jtions  (6.3.1) 
or  (6.3.2).  If  AL<  Aeff»  al  **  t°ken  as 
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Considering  the  Tower  cop  os  a  sheet  with  no  stringers,  the  effective  oreo 
may  also  be  defined,  by  the  method  outlined  in  Section  2.1.6,  at 

.  /F,  V* 

K'.i  -  (6.3.5) 


where  Per  is  the  allowable  plate  buckling  stress  given  by 


where  f(Ty2K  is  the  maximum  skin  thickness  given  by 
fmax  = 

CL 

and  Kc  is  the  plate  buckling  coefficient  (Kc  =  3.62  for  simple  supported 
edges)  or. 


KcEL  (tL\ 

Vd2j 


fc  is  the  compression  stress  on  the  lower  cop  and  is  calculated  by 

f  -Mr 
c"bwAL  ' 


(6.3.6) 


(6.3.7) 


Now  substituting  Equations  (6.3.6)  and(6.3.7)  into  Equation  (6.3.5), 
!  * 

The  maximum  stress,  fc  „  on  the  lower  cap  is  then  calculated  by 
''max  r 


(6.3.8) 


f.  =  Mr 

cmax  — — 

®w  ^eff 

where  Aeff  is  calculated  by  Equation  (6.3.8),  or 
/  rn  2  /  M-3  .  */2 


W  /  ^  Kcl 


R  ’ 


(6.3.9) 


(6.3.10) 
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Tim  value  of  fc^  is  then  compared  with  If  fCmox  -  Failow 

the  lower  cop  configuration  Is  determined  os  a  sheet  of  thickness  =  Aj/CL 
end  no  stringers.  Oft  the  other  hand,  if  fCroqx  >  Paffow  ^  fr  nccesscriy 
to  odd  stringers  to  carry  the  compression  looc!  produced  by  the  reverse  applied 
moment,  Mr.  As  an  initial  approximation  the  number  of  stringers,  n,  is 
assumed  to  be  2.  The  skin  width,  bj^  (distance  between  stringers)  is  given 

kf 


(6.T.TT) 


The  allowable  buckling  stress. 


r_  is  given  far 
a* 


where  Ki  Is  the  buckling  coefficient  for  the  skin  (K^  =  3.62  for  simple 
supported  edges).  Solving  for  the  skin  thickness,  t j_, 

1/2 

tL  =  b,  6.3.12) 


(6.3.13) 


But,  by  Equation  (6.3.3), 

Adr  -  _f?L_ 

•\vFallow 


Substituting  this  value  of  Aeff  Into  Equation  (6.3.13), 
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This  valueef  tj_  is  then  compared  with  tfie  maximum  possible  sJcin  thief: ness, 
taKB(  given  by 


t  =  AL 
Tmax  — _ 


(6.3.16) 


If  t|^  ^tfnqy,  then  n  must  be  changed  by 

«n  -  o^i+2 


(60.12V 


and  the  above  procedure  repeated  starting  with  Equation  (6.3. 111. until 
<  *max*  ^b's  procedure,  rrr  effect  adds  stringers  urrf  IT  the  cross-section 
is  able  to  sustain  the  compression  Ufcxh 

If  f|_  <tm£jx/  mus*  be  compared  with  a  minimum  skin  thickness  (.020  inch). 
If  f <  .020  t^  is  set  =  .020  and  recalculated  along  with  traax  in 
Equation  (6.3.16). 


The  allowable  compression  stress,  Fallow  may  Be  defined  by 

.  2 


Fallow-  Kcsfr  F|_ 


(6.3.18) 


where  tSfr  is  the  stringer  thickness,  bSfr  is  the  stringer  height,  and  Kc$fr  is 
the  buckling  coefficient  (KCsfF  =  *385  for  one  edge  free). 

In  terms  of  the  stringer  area,  Asfr/-  Equation  (6.3.18)  may  be  solved  for 

•sir  by 

/  2  1/4 

ttfr  =  [  FoHow  \6  1  (6.3.19) 

\KS„  EL  J 


The  stringer  area,  A^ff  may  be  defined  by 


.  _  (tmax  “  *l)  CL 

•'  ;; 


kstr  *str 


(6.3.20) 


where  tmo*  =  A[/CL  and  n  is  the  number  of  stringers.  Substituting  this 
value  of  Ajji*  into  Equation  (6,3.19), 


_  / Fallow  (fmax~fL)  CL2  \ 

*  ’  \  "2  S„Ee  ) 


(6.3.21) 
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Equation  (6.3.20)  may  be  solved  for  bj^  by 

C 


The  effective  skin  width,  beff,  (width  of  skin  that  remains  stable  under  the 
compression  load)  may  be  calculated  by 


2  (tL  beff  +  bjfr  tsfr) 


This  value  of  the  column  allowable,  Fc,  calculated  by  Equation  (6.3.24) 
must  then  be  compared  with  FaIlow.  If  Fc  <  Fa||ovv,  the  number  of  stringers 
n,  must  be  increased  as  in  Equation(6.3. 17)  and  the  above  procedure  re¬ 
peated,  starting  with  Equation  (6.3.11)  until  Fc  >  Follow* 

If  Fc  2  Fa,Jow  the  configuration  is  determined  and  the  last  values  of  Aj_, 
t|_,  t$tr,  and  n  are  taken  as  the  optimum  values. 
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LOWER  CAP 


FIGURE  6. 14 
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6.3. 1.T  DETERMINATION  OF  UPPER  CAP  AREA  AND  SKIN-STRINGB 
CONFIGURATION 

The  upper  cop  will  consist  a f  skin  and  stringers  simile*  tm  ifwt  shomt  in 
Figure  6. 15  below. 


UPPER  CAP 


FIGURE  6.15 


For  the  purpose  of  analysis  only  one  T -section  will  be  considered.  The 
temperature  of  the  upper  skin  and  stringers  are  a  function  of  the  skin 
thickness,  ts,  and  stringer  depth,  h,  respectively  as  outlined  in  Reference 
4.  A  typical  temperature-thickness  curve  is  shown  in  Figure  6.16. 

It  is  assumed  that  for  optimum  design,  the  section  will  be  critical  as  a 
column,  the  stringer  will  be  stable,  and  the  skin  will  be  subject  to  local 
buckling.  The  loading  on  the  upper  cap  is  assumed  to  be  on  applied  strain, 
eu ,  (produced  by  the  applied  moment  A^). 

For  a  given  skin  thickness,  fs,  and  stringer  spacing,  bs,  a  ratio,  R,  of 
stringer  area  to  skin  area  (R  =  tj3h/fsbj)  is  assumed  along  with  the  temperature 
of  the  stringer,  Tjj.  At  failure  the  stringer  is  ot  yield  strain,  ey^.  In  order 
to  determine  the  stringer  geomerry,  the  stringer  strain,  e^,  is  assumed  to  be 
at  the  critical  plate  buckling  strain,  ej^,  or: 
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SKIN  TEMPERATURE  VERSUS  SKIN  THICKNESS 
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eb  =  «€%=  T  (6.3.27) 

where  Kcb  is  the  buckling  coefficient  (Kc.  =  .385  for  one  edge  free),  e^i 
may  be  found  from  the  Romberg -Osgood  relationship  described  in  Section 


2.0  as 


(6.3J28) 


Solving  Equation  (6.3.27)  in  terms,  of  e^andR* 

/  T  C.  It  %  _  . 


ft*) 


(6.3.29) 


The  temperature  of  the  stringer,. Tj^/  may  now  be  found  from  figure  6.16  and 
compared  with  the  value  of  Tb  assumed;  a  new  value  of  Tb  assumed  and  the 
above  procedure  repeated  until  the  value  of  Tb  obtained  from  Figure  6.16 
agrees  with  the  value  assumed.  At  this  point  the  stringer  depth,  h,  is 
defined  for  the  assumed  value  of  R. 


The  strain  on  the  upper  skin,  ej,  is  defined  by 

«s=ds+  (cap  +  ep  +  ej)  =  e„ 

% 

where  e0  is-some  known  strain. 

The  strain  on  the  stringer,  e^,  is  calculated  by 

eb=aTb  +  (eq>  +  «p  +  eT) 


(6.3.30) 


(6.3.31) 


Solving  Equation  (6.3.30)  for  (eap  +  eD  +  ej)  and  substituting  in  Equation 
(6.3.31),  P 


eb  =  ew-KTs-Tb) 


(6.3.32) 


«s  =  eu 


(6.3.33) 


The  allowable  crippling  stress,  F  ,  as  defined  in  Section  2.1.6  is  found  by 


_  SFnE„AnCn  .  FsEstsbsCs  +  FbEb*b b 
cc  ~  ^sbs  +  Vbh 


(6.3.34) 
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where  F  is  the  stress  associated  with  the  strain  found  by  the  Ramberg-Osgood 
relationship  explained  in  Section  2.0  and  CB  is  the  effective  area  coefficient 
outlined  in  Section  Z.T.6. 

Since  the  stringer  is  a  stable  element  =  1.  On  the  other  hand,  the  skin 
element  is  subject  to  local  instabilities  and  therefore  from  Section  2.1.6, 


C*  = 


(?)  -  ?(cL) 


ft  *4  (f)' 


-  ^  es  >  «cr. 


(6.3.35) 


•  C$  =  1 


ifes  -ecrs 


The  stringer  stress,  Fb,  is  assumed  to  reach  a  cutoff  stress  where  Fb  -  F  , 
for  the  optimum  design. 


Substituting  thes  values  of  C&  and  F^  into  Equation  (6.3.34) 

t,2  CTf  (F  S)I/2  Jl+|  lJ-,/2  +  Fyb.b  h 

Ejijbj  +  Ebtbh  (6.3.36) 

The  maximum  crippling  stress  (Fcc)mox  ‘s  by  differentiating  Equation 

(6.3.36)  and  equating  the  derivative  to  zero. 


2(Est,bs  t  Ebtbh) 


<>  «y,/2  [• * 3  f»)  ]  [3 


2  (Esfsb$  +  Eb  fb  h) 


(6.3.37) 


Solving  Equation  (6.3.37)  for  Fs, 

7  I  I 

Fi  =  Fr>  3M 


(6.3.38) 
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'  Substituting  this  value  of  F  into  Equation  (6.3.36)  and  putting  t.h  in  terms 
of  R,  * 

(Fcc)max  =  r  n 

v1^  k(^)(-^rH  +F«,v  , 


(6.3.39) 


The  neutral  axis  is  located  by 


SgaAp/n,  =  LiUL 

2^*0  E«  + 


W  C.+  EuR  (2  +  t  ) 


(6.3.40) 


The  allowable  column  stress/  F  ,  is  calculated  by  Johnson’s  Parabola. 


Fc  *  (Fc=>ma*  [> 
=  (Fcc)max  + 


+  4" 2  (f/L)2^ 


a_ 

nAji_ 


(6.3.41) 


(f*cc)mcx  1-  (Q  +  C^R) 


6rr2  R 


The  elastic  applied  strain,  em,  on  the  cross-section  may  then  be  calculated 
by; 


<%>calc  *  E-SAnC„  .  Fc  (C.  +  ChR) 


-  XE„A„  E,  +  4R 

Also  eap  may  be  calculated  from  the  equilibrium  equation  by 


«ap  %,  ^rAn^r 
SEnAn 


Es  +  %R 


(6.3.42) 


(6.3.43) 


where  Fs  and  F^  are  the  stresses  associated  with  the  strains  es  and  e^  respec- 
tively,  calculated  in  Equations  (6.3.32)  and  (6.3.33). 

This  value  of  e0p  calculated  in  Equation  (6.3.43)  is  then  compared  with 
(eap)calc  found  by  Equation  (6.3.42). 

A*ap  -  eap  -  (eop)caic  (6.3.44) 

The  value  of  R  (ratio  of  stringer  area  to  skin  area)  is  then  changed  and  the 
above  procedure  is  repeated  storting  with  Equation  (6.3.29)  until  =  (eap)ca|c 
or  A e0p  =  0.  This  then  satisfies  the  conditions  of  equilibrium  ono  the 
cross-section  is  at  the  critical  column  stress  and  the  optimum  cross-section  has 
been  achieved  for  a  given  stringer  spacing.  bs,  skin  thickness,  fj,  and 
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opplled  strain*  e^. 

The  fofaf  area  of  the  upper  cap  a  then  calculated  by! 

=  &**  +  &}  =  CLt,  (1  +  R)  (6.3.45) 

6.3.}. 3  BOX  BEAM 

The  box  beam,  consisting  of  two  spor  webs,  end  upper  and  lower  caps,  h 
similar  to  that. shown  m  Figure  6v  17 below* 


FIGURE  6. 17 


The  box  beam  is  assumed  to  be  loaded  by  am  applied  moment,  M  ,  producing 
compression  in  the  upper  cap,  or  a  reverse  moment,  Mg,  producing  tension 
in  the  upper  cap,  and  a  shear  load,  Q,  on  the  webs.  It  is  assumed  that  the 
temperature  of  the  webs  and  lower  cap  are  known,  and  the  temperature  of 
the  upper  cap  is  a  function  of  the  upper  cop  thickness* 


First,  the  lower  cap  area,  Aj,  and  skin-stringer  configuration  are  determined 
by  the  method  presented  in  Section  6.3. 1.1. 

To  start,  a  value  of  the  upper  cap  stringer  spacing,  bs,  is  assumed  along  with 
values  of  the  upper  cap  skin  thickness,  t$,  and  the  upper  cop  applied  strain 


131 


NORTH  AMERICAN  AVIATION  .  INC 

COtVMtWS  tnriKON 
CtdHHVI  !♦.  »ll» 


The  web  strain,  e^,  is  calculated  by 

«w  *  « u  <VV  +  ‘?u_+ebw  (6.3.46) 

where  efe*,  is  an  assumed  bending  strain  consisting  of  both  thermal  and  applied 
bending  strains.  Tw  and  Tu  are  the  temperatures  of  the  web  and  upper  cap 
skin  (discussed  in  Section  6.3. 1.2)  respectively;  Tg  is  some  reference  temp*- 
erature  (  room  temperature);  a  u  and  a  are  the  coefficients  of  thermal 
expansion  of  the  web  and  upper  cap. 


The  lower  cap  strain,  ej_,  is  calculated  by 

CL  =  %  O’u'7^  "  +«u  +  «bL  (6.3.47) 

where  «{,.  Is  a  bending  strain  consisting  of  both  applied  and  thermal  strains 
and  Tj_  is  the  temperature  of  the  lower  cap. 

The  spar  web  strain,  ew,  h  assumed  to  be  the  strain  at  the  centroid  of  the 
web;  therefore. 


(6.3.48) 


Substituting  this  value  of  ej^  into  Equation  (6.3.46), 

«w  =  au  (Ju~Tb)  ~  °w  ^VTb)  +  eu  +  eb[/2  (6.3.49) 

and 


'1  -  *„  <VtB>  -  «l<Tl-T6>  +  to*  «bL  t6-3-50) 

A  value  of  e^  Is  assumed  (compression)  and  to  start,  let  eg^  =  0. 

From  the  strains,  eu,  ew,  and  e^,  the  stresses  associated  with  these  strains 
may  be  calculated  by  the  Ramberg-Osgood  relationship  described  in 
Section  2.0.  With  these  stresses,  f0,  fw,  and  f^,  the  corresponding  Secant 
Moduli,  Ejy,  and  Ej^,  may  be  found  by: 

Ej  =  (6.3.51) 


The  optimum  structure  Is  one  in  which  the  spar  webs  are  at  both  the  critical 
compression  buckling  stress  and  the  critical  shear  buckling  stress.  From 
this  assumption  it  is  possible  to  calculate  the  spar  web  thickness,  tw,  by  the 
method  described  In  Section  6.2.1. 
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At  tfws  point  the  upper  cap  arear  Ay,  end  the shin-stringer  configuration- 01 
determined  by  the  method  outlined  in  Section  6.3. 1.2  for  a  given  skin 
thickness,  t$,  stringer  spacing,  bjy  and  strain,  ey.  However,  if  the  value 
of  h  (stringer  height)  exceeds  some  maximum  value  (such  as  bw/2)  a  new 

_ i  _  ..f  _  it.  .! _  •  J  |f _ I _ I l.J  _  .  - _ • 


value  cf  ey  must  be  chosen  and  the  entire  procedure  started  wer  again. 


The  neutral  oasts  loco  Ted  by 

y  "  2  Eji  Aft  Yftp. 

s, 


where  Yr>f  Is  the  distance  from. some  reference  axis  to  the  centroid  of  element 
n.  In  this  cose  the  reference  axis  is  taken  to  be  an  axis  passing  through  the 
centroid  of  the  upper  skin.  Therefore, 

v  =  +  ^lA1*  rA.  3J 


(6.3.53) 


^  *w  +  A^  +  Es^  Ay 
At  this  point  the  moment  equilibrium  may  be  written  as: 
~  /TenEnAn  /n 


(6.354) 


Mu  («u  Ejy  CL  ts  Cs  +  Fjj  t]j  h)  y  -  ew  EjW  2bwtw  (  y  -  ^  ) 


+  eLEJLAL(bw-7) 


(6.355) 

I 


The  strains  ew  and  may  be  written  in  terms  of  ey  from  Equations  (6.3.49) 
and  (6.3.50).  Equation  (6.3.55)  then  becomes 


Mu  =  -<euEsuCLtsCs+Fb  fbh>  y 
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ht 

-  [«„  +  “u  <Tu  -  V  -  «W  (Tw  -  Tg)J  ESw  2bw.w  (y  -  2  ) 

+  B'+°u  <V  V  -  aLCrL-T^jeSLALo>w-7 

—  ‘  _  ^ 

+  ebL  [EsLAL  (bw-y)  -  Esw  Mw  (7"  2  )j  (6.3.5$ 

Solving  Equation  (6.3.56)  for  eb^/ 

r  ct 

eb[_=  j  Mu  +  <‘u  Es„  CL,sti  +  Fb  ^7"  *b  h)  y 

-  [«.♦  “udu-Tg)-  %  <y  -TT ) 

* 

“[eu  +  ou(T0-TB)-  qL(Tl-Tb)]  [EsLAL(bw-7) 


+  Eswtwbw(y”  2 


(6.3.57) 


Trie  value  of  eb^  calculated  by  Equation  (6.3.57)  is  then  compared  with  the 
value  of  eb^  assumed  in  Equations  (6.3.49)  and  (6.3.50).  A  new  value  of 
ebb  ‘s  assumed  and  the  steps  outlined  in  Equations  (6.3.49)  through  (6.3.57) 
are  repeated  until  the  value  of  ebb  calculated  by  Equation  (6.3.5 7)  is  com¬ 
patible  with  the  value  c!  ebb  assumed. 


Next,  the  axial  strain  equilibrium  equation  may  be  written. 
XenEnAn  =  0 


(6.3.5^ 


eu  ^Sy  *s  +  bs  *b  b  +  2  ew  Nvbw  +  el  \  ® 

where  is  the  stress  associated  with  the  strain  on  the  stringer  calculated  in 
the  upper  cap  in  Section  6. 3. 1.2.  Solving  Equation  (6.3.58)  for  ey/ 


“  ^b  £L_  tjj  h  -  2  ew  E^  twbw  -  e^  Ej^ 

K _ __________ 

EsyCsClt, 


(6.3.59) 
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The  value  of  eu  calculaled  by  equation  (6.3.  5^  Is  then  compared  with  the 
value  of  eu  assumed  in  Equation  (6.3.49}:  A  new  vafue  of  eu  is  assumed  and 
the  above  procedure,  starting  with  Equation  (6.3.49),  is  repeated  until  the 
value  of  e0  calculated  by  Equation  (6.3.59)  is  compatible  with  the  assumed 
value  of  e0. 

The  structure  is  now  in  complete  equilibrium  since  the  conditions  imposed  hy 
Equations  (6.3.57)  and  (6.3.59)  have  been  satisfied.  . 

Next,  a  new  value  of  upper  skin  thickness,  ts,  is  assumed  and  the  entire 
procedure  outlined  above  is  repeated  until  the  necessary  equilibrium  con¬ 
ditions  are  satisfied.  After  several  values  of  the  upper  skin  thickness  have 
been  assumed,  a  plot  may  be  made  of  the  total  area,  Ay  (Ay  =  A  0+  2bwtw 
+  Ai),  versus  ts  for  a  given  stringer  spacing,  bs,  os  slrown  in  Figure  6.18 
(a).  This  curve  exhibits  a  minimum  value  of  Ay  at  some  value  of  t$.  This 
point  is  the  optimum  design  for  a  given  stringer  spacing  bs.  The  minimum 
point  may  be  calculated  by  determining  three  point,  and  assuming  a  parabola 
through  them  as  discussed  in  SecSion  6.4. 

A  new  value  of  the  stringer  spacing,  bs,  is  assumed  and  the  dbove  proceckire 
repeated  for  several  values  of  ts  to  determine  another  curve  of  Ay  versus  t$. 

A  typical  graph  of  Ay  versus  t  for  several  different  values  of  bs  is  shov.n  in 
Figure  6. 18  (b).  A  plot  Is  then  made  of  (Ay)mjn  versus  bs  as  shown  in  Figure 
6.18(c). 

The  stringer  spacing,  bs,  at  the  minimum  of  the  curve  shown  in  Figure  6.18 
(c)  is  the  optimum  stringer  spacing  and  the  optimum  structure  has  been  defined. 
However,  it  is  possible  that  the  upper  cap  may  fail  in  tension  when  the 
reverse  moment,  Mr,  is  applied.  Thus  the  upper  cap  area  must  be  tested  to 
see  if  it  is  capable  of  sustaining  such  a  tension  load  and,  if  necessary, 
adjust  the  upper  cap  area  and  skin  stringer  configuration  to  sustain  the  tension 
load  produced  by  Mr.  The  minimum  upper  cap  area,  Amin,  required  to 
sustain  a  tension  load  produced  by  the  reverse  moment,  Mj^,  .is  calculated  by 


(6.3.60) 

where  Ft„  and  FtWL  are  the  ultimate  tension  stresses  of  the  uppc.'  skin  and 
stringers  respectively.  If  the  minimum  area  shown  in  Figure  6. 18(d)  is  less 
than  Amjn  calculated  by  Equation  (6.3.60)  a  new  bs  must  be  chosen  such  that 


%  (Es  »s  bs  +  Eb  tb  h) 


F*us  Es  *S  bs  + 
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£  roin  -  end  the  enlire  procedure  ii  repeated  one  «oe». 

time  for  this  new  value  of  bs.  This,  then  is  this  final  optimum  design  of  the 
box  beam.  It  should  be  noted,  however,  that  this  is  an  optimum  design  for 
a  given  rib  spacing,  L,  chord  length,  CL,  box  depth,  by,,  applied  primary 
moment.  My,  reverse  moment,  shear  load,  Q,  and  temperature  distribu¬ 
tion  through  the  spar  webs  and  lower  cap. 

The  IBM  program  to  perform  this  optimization  Is  presented  in  Appendix  D. 
6.3.2  BOX  BEAM  WITH  CORRUGATED  WEBS' 


It  was  pointed  out  In  Section  6.2.3  that  the  corrugated- spar  web  design 
should  be  considered  when  thermal  stresses  are  present  since  the  corrugated 
spar  web  will  not  restrain  the  hot  upper  cop  from  expanding  and  the  thermal 
stresses  induced  In  the  structure  ore  much  lower,  thus  permitting  a  lighter 
weight  design.  However,  the  comigated  spar  webs  will  allow  greater  de¬ 
flections  and  a  considerable  weight  penalty  must  be  paid  in  the  case  of  a 
stiffness  design. 


The  load,  P,  carried  by  the  upper  and  lower  caps  can  be  determined  directly 
for  the  corrugated  spar  web  case.  P  may  be  calculated  by: 


(6.3.61) 


The  upper  cap  area.  Ay,  and  skin-stringer  configuration  is  determined  by  the 
method  described  in  Section  6.3.1. 2  after  a  value  is  assumed  for  the  upper 
cap  applied  strain  eu.  P  may  then  be  calculated  by: 

P  =  fstsCL+fb^|:  tyk  (6.3.62) 

where  fs  and  are  the  stresses  associated  with  the  strains  on  the  upper  skin 
and  stringers  respectively.  The  load  calculated  by  Equation  (6.3.62)  is  then 
compared  with  the  actual  load  calculated  in  Equation  (6.3.61)  and  a  new  value 
of  eu  is  chosen  and  the  procedure  repeated  until  the  calculated  load  is  com¬ 
patible  with  the  actual  load. 


The  lower  cap  area,  Aj,  and  skin-stringer  configuration  is  determined  by  the 
method  outlined  in  Section  6.3. 1.1. 


The  design  of  the  corrugated  web  is  then  determined  by  the  procedure  described 
in  Section  6.2.1. 
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The  design  procedure  to  this  point  has  been  for  a  given  upper  skin  thickness,, 
tj,  and  stringer  spacing,,  As.  in  the.  straight  weixcaser  two  more  veives-ef 

fs  are  chosen  for  a  given  bs.  The  orecr  of  each  configuration  is  computed  by: 

Aj  -  Ay  +  2  A w  +  A^  (6.3.63) 

The  parabola  method  of  Section  6.4  is  again  used  to  determine  the  minimum 
area  configuration  as  a  function  of  the  skin  thickness,  Cv  This  procedure  uc 
repeated  for  two  more  values  of  b$  to  determine  the  minimum  area  configuration 
with  respect  to  ts  and  bj. 

The  IBM  program  for  both  the  straight  ond  corrugated  spar  webs  Is  presented 
in  Appendix  D. 
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6.4  PARABOLA  CURVE  FITTING  METHOD 

The  porobofa  curve  fitting  method  rirused  in  this  study  for  tfrepiiewiy 
purpose  of  describing  the  optimization  curves  generated  by  the  high  speed 
computer  data.  The  general  form  of  the  equation  used  in  this  method  Is: 

y=Ax2  +  Bx  +  C  (6.4.1) 

To  compute  the  exact  equaffon  for  the  optimization  curves,  a  minimum  of 
three  paints  are  required.  These  three  points  ccn  be  described  as  xjy^, 
*272'  cn£^  x3>«3.  Using  these  points  in  the  general  equation,  the  following 
three  equations  can  be  written: 

r 

yj  =  Axj  +  Bxj  +  C 

y2  =  Ax*  +  Bx2  +  C  (6.4.2) 

2 

y3  =  Ax3  +  Bx3  +C 

Solving  these  three  equations  simultaneously  the  following  values  for  the 


constants  A,  B,  and  C  are  obtained: 

A  = 

(y?  -  yq) 
(*•?  -  *a) 

(*2  +  *3^ 

-  (dzyA 

(x1  -  x9) 

-  •  (xj  +x2) 

(6.4.3) 

B  = 

yi  -  y2 
xi  -  x2 

A 

X1  +  x2 

(6.4.4) 

C  = 

y3  - A  x3  + 

(yi-y2)  ‘  Ax3 

(X1  -  x^  x,  -  x2 

(6.4.5) 

Substituting  these  constants  into  the  general  equation,  an  equation  through 
the  three  chosen  points  is  obtained.  In  this  study  the  primary  objective  is 
to  estimate  the  minimum  value  of  the  parameter  y  with  respect  to  x.  The 
point  on  the  curve  at  which  this  minimum  occurs  has  a  slope  equal  to  zero. 
From  calculus  it  is  known  that  the  first  derivative  of  an  equation  gives  the 
slope  at  any  point  on  the  curve.  The  derivative  of  the  general  equation  is: 

=  2  Ax  +  B  (6.4.6) 

If  the  slope  dy/dx  is  set  equal  to  zero  the  equation  can  be  solved  for  the 
parameter  x.' 
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0  -  2  A  x  +  B  (6.4.7) 

~2A 


This  value  of  x  is  then  the  point  at  which  the  minimum  or  maximum  value  of 
y  occurs.  Certain  controls  have  been  placed  in  the  computer  program  to 
insure  that  it  is  a  minimum  which  is  computed  and  that  the  value-iso 
relatively  accurate  vaioev 
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6.5  CONCLUSIONS 

hr  summary,  optimization  methods  have  Seen  developed  for  four  types  d  box 
beams;  honeycomb  box  beam  with  straight  or  corrugated  webs  and  Tntegrdf 
sking-sfringer  box  beam  with  straight  or  corrugated  webs.  One  of  the 
principle  features  of  these  optimization  methods  is  their  ability  to  handle 
variable  temperature  gradients.  However,  as  with  all  programs,  certain 
restrictions  and  limitations  must  be  imposed.  The  following  paragraphs  describe 
some  of  these  restrictions  end  limitations  along  with  soma  suggested  appli¬ 
cations  of  the  programs. 

U  m  ist  be  pointed  out  the  optimization  program  presented  tn  Section  6  is  f or 
a  variable  mold  line  since  the  bole  beam  depth,  is  taken  as  the  distance 
between  the  centroids  of  the  upper  and  lower  cops.  Thus  the  mold  line  will 
vary,  depending  on  the  skin  thickness  and  stringer  depth.  However,  this 
variation  in  mold  line  will  be  small  and  a  good  approximation  may  be  made 
of  the  optimum  design. 

The  optimization  methods  presented  here  are  for  a  fixed  rib  spacing,  L. 
However,  the  designer  may  easily  optimize  for  the  rib  spacing  by  calculating 
the  optimum  design  for  several  different  rib  spacings  and  including  an 
additional  term  to  account  for  the  weight  of  the  rib  as  a  function  of  the  rib 
spacing. 

The  optimum  design  at  elevated  temperatures  may  not  be  optimum  at  room 
temperature  and  may,  in  tact,  not  be  able  to  sustain  the  applied  load  at 
room  temperature.  Therefore,  the  programs  presented  here.must  be  run  for 
room  temperature  to  see  If  the  box  beam  is  indeed  critical  at  room  temperature. 
A  plot  may  be  made  of  the  room  temperature  optimum  curve  as  shown  in 
Figure  6. 19.  The  intersection  of  these  two  curves  is  then  the  optimum  design 
point  when  both  room  temperature  and  elevated  feirperature  are  considered, 
except  when  the  minimum  of  one  of  the  curves  lies  inside  the  envelope 
described  by  the  two  curves. 


SUPERPOSITION  OF  ROOM  AND  ELEVATED  TEMPERATURE 
OPTIMIZATION  CURVES 

Figure  6. 19 
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This  principle  of  superposition  of  optimum  curves  may  be  extended  to  other 
para  maters  proceed  by  different  flight  conditions  and  o  design  envelope 
established  to  determine  on  optimum  structure  that  will  satisfy  off  of  the 
conditions  imposed  upon  it.  A  typical  example  is  shown  in  Figure  6.20 
below. 


Condition  1 


'  Design  Variable 


SUPERPOSITION  OF  MULTI-CONDITION  OPTIMIZATION 
CURVES 


FIGURE  6.20 


